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SECTION  I 


INTRODUCTION 

NELSIM  (Nonelectrical  Languages  Simulation  Module)  is  a digital  com- 
puter program  written  in  FORTRAN  IV.  The  program  translates  non-electrical 
system  descriptions  into  electrical  analogs  and  differential  equations 
acceptable  for  execution  by  the  SCEPTRE  (System  for  Circuit  Evaluation  and 
Prediction  of  Transient  Radiation  Effects)  simulation  program.  The  program 
thus  provides  an  interface  between  system  level  problems  and  electrical 
simulation  programs.  The  module  allows  engineers  not  familiar  with 
disciplines  generally  associated  with  circuit  simulation  programs  to  solve 
system  problems  without  manually  deriving  the  electrical  analogs  and 
equations  necessary  to  run  such  programs.  The  input  language  allows  input 
of  mechanical , thermal,  electro-mechanical  and  electro-optical  problems  in 
terms  of  differential  equations,  transfer  functions,  or  elements  and  units 
familiar  to  those  particular  disciplines. 

The  program  features  a f^ee-field,  user-oriented  input  language 
modeled  after  the  SCEPTRE  input  language.  The  output  is  also  in  terms 
of  the  SCEPTRE  language  and  allows  direct  execution  by  SCEPTRE.  The  output 
of  NELSIM  can  be  punched  cards  if  requested  by  the  user.  Another  option 
is  scaling  which  allows  the  user  to  obtain  an  electrical  analog  with  the 
units  scaled  to  the  actual  response  magnitudes  expected. 

1.1  NELSIM  Input/Output  Options 

The  program  allows  the  three  types  of  input  to  describe  system  problems 
illustrated  in  Figure  1. 

The  transfer  function  input  allows  users  to  enter  the  description  of 
a system  or  part  of  a system  in  terms  of  signal  flowgraphs  or  block 
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DIFFERENTIAL  EQUATIONS 


TRANSFER  FUNCTIONS 


Figure  1.  Input/Output  of  NELSIM 
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diagrams.  The  language  for  this  option  is  nearly  identical  to  the  input 
language  of  the  program  SYNAP  (Symbolic  Network  Analysis  Program)  developed 
for  The  Air  Force  Weapons  Laboratory  (AFWL)  in  1972  and  consists  of  descrip- 
tions of  S-Domain  transfer  functions  connected  bwtween  user  derived  node 
names . 

The  differential  equations  option  allows  the  user  to  enter  problem 

X L. 

definitions  in  terms  of  sets  of  n . order  differential  equations  algebrai- 
cally coupled  together. 

The  user  may  also  enter  problems  in  terms  of  functional  lumped  elements 
and  nodal  connections.  The  program  accepts  elements  of  mechanical , thermal, 
electro-mechanical  and  electro-optical  systems. 

The  output  of  the  program  consists  of  a set  of  first  order  differential 
equations  for  all  differential  equation  entries  and  all  transfer  function 
entries.  It  also  consists  of  electrical  lumped  element  equivalent  networks 
for  all  entires  containing  system  elements.  The  user  may  use  any  option  or 
combination  of  options  to  define  the  overall  system  to  be  analyzed.  The 
user  may  have  the  output  of  NELSIM  punched  on  cards  and  thus  generate  an 
input  deck  to  the  SCEPTRE  program  directly.  The  user  may  also  call  on  the 
scaling  algorithm  built  into  NELSIM  and  obtain  a listing  of  the  analog 
generated  with  constant  element  values  scaled  to  improve  numerical  accuracy 
upon  execution. 

1.2  NELSIM  Input  Data  Organization 

The  data  organization  and  input  language  of  NELSIM  have  been  kept 
identical  to  that  of  SCEPTRE  wherever  possible  and  consistent  with  it  in 
all  other  cases.  Figure  2 illustrates  the  data  organization.  As  shown, 
the  data  are  entered  following  heading  and  subheading  cards.  All  problem 
descriptions  entered  in  terms  of  lumped  elements  are  entered  under  the 
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MECHANICAL  

r ELEMENTS 


Functional 
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MAIN  
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DEFINED  PARAMETERS 
FUNCTIONS 

INITIAL  CONDITIONS 
L OUTPUTS 

THERMAL  


ELEMENTS 


OUTPUTS 

ELECTRO  MECHANICAL  

ELEMENTS  (Mechanical) 
ELEMENTS  (Electrical) 
MODELS 

DEFINED  PARAMETERS 
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INITIAL  CONDITIONS 
OUTPUTS 

ELECTRO  OPTICAL  
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OUTPUTS 

DIFFERENTIAL  EQUATIONS 
COMMENTS 
EQUATIONS 
DEFINED  PARAMETERS 
FUNCTIONS 
OUTPUTS 

TRANSFER  FUNCTIONS 
COMMENTS 
BRANCHES 

DEFINED  PARAMETERS 


OUTPUTS 
RUN  CONTROLS 
RERUN  DESCRIPTION  (M) 
ELEMENTS 

DEFINED  PARAMETERS 
RUN  CONTROLS 

■END 


Main 

Subhead!  ngs 


Figure  2.  NELSIM  Input  Data  Organization 
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ANALOG  DESCRIPTION  heading.  All  sets  of  differential  equations  are  entered 
under  the  DIFFERENTIAL  EQUATIONS  heading  and  all  transfer  function  descrip- 
tions are  entered  under  the  TRANSFER  FUNCTIONS  heading.  The  RUN  CONTROLS, 
RERUN  DESCRIPTION  and  END  headings  control  the  execution  of  the  run, 
parameter  value  changes  and  program  termination,  respectively.  The  function- 
al subheadings  shown  in  Figure  2 identify  the  type  of  system  being  entered 
if  lumped  networks  are  entered.  The  functional  elements  serve  to  identify 
the  remainder  of  variable  values  such  as  element  values  and  description  of 
equations  used.  The  main  headings  must  appear  in  the  order  shown  in  the 
figure.  Main  headings  as  well  as  subheadings  however  need  not  be  included 
if  not  utilized.  For  example,  if  a problem  is  described  strictly  in  terms 
of  differential  equations,  the  ANALOG  DESCRIPTION  headings,  TRANSFER  FUNCTION 
heading  and  all  related  subheadings  may  be  deleted  from  the  input.  Sections 
II  through  IV  elaborate  the  details  of  the  input  language  as  applicable  to 
each  input  type. 

1.3  NELSIM  Usage  Highlights 

This  paragraph  provides  the  user  a quick  introduction  to  the  input 
language  by  summarizing  the  input  through  three  sample  problem  inputs.  As 
discussed  above  three  types  of  input  are  allowed  in  NELSIM.  The  input  data 
organization  for  a differential  equations  problem,  the  input  organization 
for  a transfer  function  problem,  and  system  el ement  entry  for  a mechanical 
problem  are  shown  in  Figure  3. 

The  input  data  consists  of  descriptive  statements  constructed  from 
user-derived  component  names,  node  names  and  value  specifications.  The 
information  is  placed  on  cards  separated  by  special  characters  such  as  a 
comma,  dash,  parenthesis,  or  equal  sign.  In  general,  several  complete 
statements  can  be  punched  on  a card,  separated  only  by  a comma.  The 
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Equations 
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NELSIM  J 
Input  S 


r 


K!  ( V2“V-i ) + M2 


+ k2(v2-v3)  = 0 


d2V. 


K2(V3"V2>  + M3  Hr  + K3V3  = 0 


Parampter 

values 


^ = cos  2t 


DIFFERENTIAL  EOUATIONS 
EQUATIONS 

PK1*(PV2-PV1)+PM2*D2PV2+PB*(DIPV2-D1PV3)+PK2*(PV2-PV3)=0 

PB*(D1PV3-D1PV2)+PK2*(PV3-PV2)+PM3*(D2PV3)+PK3*PV3=0 

DEFINED  PARAMETERS 

PK1=1 

PK2=2 

PM2=.2 

PM3= . 3 

PB=1 

PVI=Q1 (TIME) 

FUNCTIONS 

Q1  (T)  = (C0S(2.*T)) 

OUTPUTS 

PV2(V2) ,PV3(V3) .PLOT 
RUN  CONTROLS 
STOP  TIME=100 

INTEGRATION  ROUTINE=IMPLICIT 
END 


a.  Differential  Equations 


Figure  3.  NELSIM  Inputs 
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(input) 


(output) 


TRANSFER  FUNCTIONS 
BRANCHES 
1-2(1. ,0) 

2- 3  3.  ,1,  1 . .0/  .2,3,  1.2,2,  l.,l) 

3- 2( -1 . ,0) 

DEFINED  PARAMTERS 

Input  at ^Pl  =T1 

Node  1 FUNCTIONS 

T1=0,0,  0,1,  10,1 
OUTPUTS 

Pl,P3,PL0T-c Output  requests  at 

RUN  CONTROLS  Nodes  1 and  3 

STOP  TIME=10 

INTEGRATION  ROUTINE  = IMPLICIT 
END 


b.  Transfer  Function 
Figure  3.  NELSIM  Inputs  (Cont'd) 
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ANALOG  DESCRIPTION 
MECHANICAL 
ELEMENTS 
VI ,G-3=Q1 (TIME) 

K1 ,G2-3=1 

K2,G1-2=P2 

K3,G-1=3 

M2,2-G2=.2 

M3,1-G1=.3 

BI  ,G1 -2=1 

FUNCTIONS 

Q1  ( T ) = ( COS ( 2 - *T ) ) 

DEFINED  PARAMETERS 

P2=2 

OUTPUTS 

VM2( V2) ,VM3(V3) ,PLOT 
RUN  CONTROLS 

INTEGRATION  ROUTINE=IMPLICIT 

STOP  TIME=10 

END 


c.  Mechanical  Analog 


Figure  3.  NELSIM  Inputs  (Cont'd) 
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information  is  all  placed  under  major  heading  and  subheading  such  as 
TRANSFER  FUNCTIONS  and  DEFINED  PARAMETERS. 

1.4  Manual  Organization 


This  section  has  briefly  described  the  capabilities  of  the  NELSIM 
program  and  its  general  usage,  and  has  given  some  examples  of  the  input 
language.  The  remainder  of  the  manual  expands  on  the  topics  discussed 
above.  Section  II  discusses  the  input  language  applicable  to  non-electri- 
cal system  topologies.  Sections  III  and  IV  describe  the  input  language 
applicable  to  differential  equations  and  transfer  functions,  respectively. 
The  final  section,  Section  V,  contains  several  example  problems  varying  in 
complexity  and  mode  of  analysis.  These  examples  have  been  included  to 
help  the  reader  comprehend  the  wide  usage  range  of  NELSIM. 
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SECTION  II 

ANALOG  SYSTEM  INPUT  LANGUAGE 

As  discussed  in  Section  I,  NELSIM  allows  the  input  of  nonelectrical 
systems  in  terms  of  networks  made  up  of  lumped  elements  normally  utilized 
within  a particular  discipline.  The  disciplines  allowed  are  mechanical, 
thermal,  electromechanical  and  electro-optical . The  general  input  language 
for  the  analog  description  is  discussed  below  followed  by  the  documentation 
of  each  Individual  discipline. 

2.1  GENERAL  ANALOG  SYSTEM  INPUT  LANGUAGE 

All  problems  to  be  entered  utilizing  lumped  elements  must  be  entered 
using  a nodal  topology  description.  The  general  features  of  the  network 
input  language  are  explained  using  thermal  quantities  and  elements.  The 
entries  covered  below  are  applicable  to  all  the  systems  disciplines 
allowed.  Section  2.2  through  2.5  describe  the  input  of  mechanical,  thermal, 
electromechanical  and  electro-optical  systems,  respectively. 

2.1.1  Network  Preparation 

For  a portion  of  a system  described  in  terms  of  thermal  elements, 
the  network  topology  is  accepted  by  NELSIM  in  terms  of  the  standard  elements 
shown  in  Table  I.  Each  element  may  be  linear  or  nonlinear.  The  six  steps 
to  be  followed  In  preparing  a network  for  a NELSIM  run  are  summarized 
below: 

• Draw  an  equivalent  network  comprised  of  thermal 
resistors,  thermal  capacitors,  and  temperature 
and  heat  flow  sources 

• Assign  a name  or  number  to  all  nodes  in  the  circuit 
t Give  a name  to  each  circuit  element  (beginning  with 

R,  C,  T,  Q) 
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• Assume  arbitrary  heat  flow  directions  in  each  passive  circuit 
element 

• Indicate  the  direction  of  positive  heat  flow  In  each 
temperature  and  heat  source 

• Choose  and  record  circuit  values  in  a consistent  set  of 
parameter  units. 

A node  is  defined  as  the  point  of  common  potential  (temperature)  at  the 
junction  created  by  the  connection  of  two  or  more  network  elements.  Each 
node  is  given  an  arbitrary  name  consisting  of  six  alphanumeric  characters 
or  less.  Element  names  are  chosen  arbitrarily  but  must  begin  with  the 
letter  specified  in  Table  I and  must  not  contain  more  than  5 characters. 


TABLE  I.  STANDARD  THERMAL  ELEMENTS 


ELEMENT 

SCHEMATIC  REPRESENTATION 

PROGRAM  NOTATION 

Resistor 

If 

R 

Thermal 

Capacitor 

V\A/ 

C 

Temperature 

Source 

T 

V_7 

Heat  Flow 
Source 


2.1.2  Element  Input  Description  and  Assumed  Polarity 

In  preparing  the  circuit  a heat  flow  direction  and  temperature 
difference  are  assumed  for  each  element  as  shown  in  Figure  4. 


FLOW  DIRECTION 

Figure  4.  Assumed  Element  Heat  Flow  Direction 


The  elements  are  input  into  the  program  following  the  subheading  card 
ELEMENTS  using  the  format: 


where  the  heat  flow  direction  is  assumed  to  be  from  Node  1 to  Node  2 and 
the  parame.ers  are  defined  as  follows: 

• ELEMENT  - The  name  of  the  element  beginning  with 

the  letter  assigned  in  Table  I and 
followed  by  4 alphanumeric  characters 
or  less  serving  as  identifiers,  i.e., 

(R12,  TIN). 

e NODE  1 - Initial  node  connected  to  the  element 
and  denoted  by  up  to  six  alphanumeric 
characters.  (From  node) 

• NODE  2 - Final  node  connected  to  the  element  and 

denoted  by  up  to  six  alphanumeric 
characters.  (To  node) 
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• VALUE  - Value  of  the  element.  The  form  of  VALUE 
may  be  expressed  in  six  different  formats 
to  allow  for  linear  and  non-linear  values 
of  elements.  The  six  formats  available 
are  discussed  in  Section  2.1.4. 

The  Information  may  be  placed  anywhere  on  the  card  between  columns  1 and 
72  with  the  delimiters  (cornna,  dash,  equal  sign)  placed  in  the  order 
shown.  More  than  one  element  may  be  Included  on  any  card  if  separated 
by  coronas. 


The  heat  in  a source  Is  assumed  to  flow 
from  the  negative  to  positive  pole  and 
VALUE  Is  positive  If  the  nodal  connection 
Is  entered  with  the  "FROM"  node  denoting 
the  negative  terminal.  The  source  in 
Figure  5 for  example  can  be  entered  in 
the  two  ways  shown  below. 


/ Tl,  1-2  = 100 


/ Tl,  2-1  = -100 


Figure  5.  Temperature  Source  Polarity  and  Assumed  Heat  Flow 


Example 

The  element  entry  Into  the  program  for  the  simple  schematic  shown 


Is  given  below: 


/ TIN.  GND-1  = VALUE 
/ Rl.  1-A1  ■=  VALUE 

/ Cl.  Al-GND  - VALUE 
/ RA.  A1-2  = VALUE 

/ QOUT,  2-GND  = VALUE 


0 Rl  A1  RA  @ 
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An  equally  valid  entry  Is: 


/TIN,  GND-1  = VALUE,  R1 , 1-A1  = VALUE 
/H,  Al-GND  = VALUE  ~ 

/ RA,  Al-2  = VALUE,  QOUT,  2-GND  = VALUE 

2.1.3  Element  Temperatures,  Heat  Flows 

The  temperature  or  heat  flow  associated  with  any  element  is  referred 
to  by  placing  a T or  an  Q before  the  element  name.  TR1  and  QR1  therefore 
refer  to  the  temperature  and  heat  flow  associated  with  resistor  R1 . 

2.1.4  Element  and  Circuit  Value  Formats 

An  element  value  or  circuit  quantity  may  be  expressed  in  the  six 
formats  listed  below.  Each  is  explained  in  detail  in  the  following  para- 
graphs. 

• Defined  Parameter 

• Constant 

• Table 

• Equation 

t FORTRAN  Function  Subroutine 

• Expressions 

2. 1.4.1  Quantities  Expressed  as  Defined  Parameters 

A circuit  element  or  circuit  quantity  (temperature,  heat  flow,  deri- 
vative) may  be  expressed  In  terms  of  a DEFINED  PARAMETER.  The  format 
for  an  element  input  Is 

ELEMENT,  NODE  1-NODE  2 = PX 

Where  X Is  the  name  assigned  to  the  parameter  and  may  contain  up  to  five 
alphanumeric  characters.  The  prefix,  P,  must  always  proceed  the  name  X 
as  follows: 
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/Rl,  1-2  = P12 
/RSHNT , NOUT-GND  = POUT 

The  actual  values  of  defined  parameters  used  are  defined  under  the  sub- 
heading DEFINED  PARAMETERS.  This  Input  and  its  related  formats  are  dis- 
cussed in  paragraph  2.1.5.  A defined  parameter  may  be  expressed  in  terms 
of  a constant,  table,  equation,  function  subroutine  or  another  defined 
parameter. 

2. 1.4. 2 Quantities  Expressed  as  Constant 

An  element  value,  circuit  quantity  or  defined  parameter  may  be 
expressed  as  a numerical  constant  written  in  either  integer  or  decimal  form 
with  or  without  an  exponent.  Up  to  13  characters  may  be  used  to  represent 
the  constant.  Table  II  summarizes  the  constant  forms  acceptable  to  the 
program.  Sample  element  description  cards  using  the  constant  notation  are 
illustrated  below. 

Table  II.  Acceptable  Constant  Formats 

/Rl,  1-2  = 10 
/R2 , 3-4  = l.E-3 
/C6,  GND-1  = E-12 


2. 1.4. 3 Quantities  Expressed  as  TABLES 

A variable  circuit  quantity  or  defined  parameter  may  be  expressed 
In  tabular  form.  The  table  format  is  specified  as  follows  for  an  element: 

/ELEMENT,  NODE  1 - NODE  2 = TABLE  X (VARIABLE) 


9*' 

Format 

Example 

Value 

+ XX 

10 

10. 

+ XX. 

-10. 

-10. 

+ X.  X 

-10.  2 

-10.  2 

+ X.  XE+Y 

+ 10.  2E-6 

10.  2 x 10*6 

+ XXE+Y 

-10E-4 

-10.  X10"4 
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/ELEMENT,  NODE  1 - NODE  2 = TX  (VARIABLE) 


where  X is  the  table  name  containing  up  to  five  alphanumeric  characters 
and  the  variable  in  parenthesis  is  the  independent  variable.  The  inde- 
pendent variable  may  be  any  circuit  variable  such  as  voltage,  current, 
elements  and  defined  parameters. 

Source  Q1  is  a function  of  the  temperature  across  capacitor  Cl  as 
shown  in  the  table 


Ql 

TCI 

0. 

0 

1. 

1 

1. 

2 

1.5 

2 

1.5 

3 

0. 

3 

0. 

4 

Q1  could  be  entered  as 


/Ql,  1-2  = TABLE  KTC1!) 


or 

/Ql,  1-2  = Tl(TCl) 

The  tabular  data  shown  In  the  table  is  entered  as  Input  pairs  under  the 
subheading  FUNCTIONS  and  Its  format  Is  discussed  In  Section  2. 1.4. 6.  If 
the  parentheses  and  variable  are  left  off  the  table  description,  TIME  is 
assumed  to  be  the  Independent  variable.  Sample  table  entries  follow. 


16 


/R12A,  A1-A2  = TABLE  T1 ( QR1 


/RA 14,  8-12  = T3(TC3] 


/TIN,  GND-I  = TTIN  •* — Time  assumed  as 

independent  variable 

/Q1 , 2-3  = T14(P3) 

The  occasion  often  arises  that  the  same  table  data  is  utilized  by  more 
than  one  element  or  parameter.  A common  situation  occurs  when  two  current 
generators  in  the  network  are  defined  by  the  same  tabular  data,  differing 
only  in  the  independent  variable.  In  this  case  the  user  may  define  both 
generators  under  ELEMENT  as 

/Ql,  1-2  = TABLE  1 (TCI ) 

/Q2,  7-8  = TABLE  1 (TC2) 

2. 1.4. 4 Quantities  Expressed  as  Expressions 

A circuit  quality,  element  or  defined  parameter  may  be  expressed 
in  terms  of  an  expression.  The  element  input  format  for  this  form  is: 


LEMENT,  N0DE1 -N0DE2  = XN  (NUMERICAL  EXPRESSION 

where  N is  the  expression  denotation  containing  up  to  five  alphanumeric 

characters.  An  example  is  given  below  for  source  T1 
1 

9 Al.G-l  = XT  (10.*C0S(2.*3.14*PF) 


T1  = 10. COS  2irf  where  PF  = f 


2. 1.4. 5 Quantities  Expressed  as  EQUATIONS 


A circuit  quantity,  element  or  defined  parameter  may  be  expressed 
in  terms  of  an  equation.  The  element  input  format  for  this  form  is: 
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/ELEMENT,  NODE  1-NODE  2 = EQUATION  X(a] , a2>  . . . ap  ) 

or 

/ELEMENT,  NODE  1-NODE  2 = QX(ar  a?,  . . . a ) 

where  X is  the  name  containing  up  to  5 alphanumeric  characters  given  to 
the  equation.  The  variable  a-j  through  ap  are  described  in  the  equation 
and  may  be  any  circuit  quantity,  number,  defined  parameter  or  table. 

Example 

Resistor  R3  between  nodes  1 and  2 is  defined  as: 

R3  = f (QR3,  TCI,  P4) 

where  f (QR3,  TCI,  P4)  = (QR3  - TC1)*P4. 

The  element  input  entry  could  be  as  follows: 

/R3,  1-2  = EQUATION  1A(QR3,  TCI,  P4) 

OR 

/L3,  1-2  = Q1A(QR3,  TCI,  P4) 

The  actual  mathematical  expression  f(QR3,  TCI,  P4)  = (QR3  - TCI )*P4  must 
be  defined  under  the  subheading  card  FUNCTIONS.  The  details  of  this  input 
are  discussed  in  Section  2.9.  More  sample  inputs  are  provided  below. 

/C2,  6-1  = EQUATIONS  2( . 5,  120.,  TABLE  4(TC2)) 

/m/,  1-12  = Q3(TR3,  TR2VT2.,  Tl(0&2)) 

NOTE:  DECIMAL  POINTS  MUST  BE  SUPPLIED  FOR  ALL  CONSTANT 

QUANTITIES  ENTERED  IN  THE  ARGUMENT  LIST  at  THROUGH 

a . 

-n 
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2. 1.4. 6 Quantities  Expressed  as  FORTRAN  FUNCTION  SUBROUTINES 

The  user  knowledgeable  of  FORTRAN  may  enter  circuit  elements, 
quantities  and  defined  parameters  as  separate  FORTRAN  function  subroutines. 
This  feature  is  extremely  useful  when  a quantity  is  too  complex  to  be 
expressed  as  an  equation.  The  element  input  format  for  this  form  is 

/ELEMENT  NAME,  NODE  1 -N0DE_2  = FUNCTION  X(a]t  a?,  . . . , aj 

or  simply 

/ELEMENT  NAME,  NODE  1 NODE  2 = FX(a1 , a?,  . . . , a ) 

where  X or  FX  is  the  name  of  the  subroutine  and  contains  up  to  five  alpha- 
numeric characters  and  a-j  through  an  are  the  variables  utilized  in  the 
function  subroutine.  These  variables  may  be  any  circuit  quantity,  element 
or  defined  parameter.  The  subprogram  name  must  be  unique  and  begin  with 
the  letter  F if  the  short  form  above  is  used. 

Example 

Capacitor  Cl  between  nodes  1 and  2 is  defined  as  follows: 

Cl  = f (TIN,  TC3)  for  QC1  IE-3 

Cl  = f (TOUT , TC3)  for  QC1  IE-3 

where  f (TIN,  QC3)  = 

ctTr\  iit  Tr**\  _ T0UT-TC3 
and  f(T0UT,  TC3)  = — — 

Using  the  subroutine  input  option  the  entry  is  as  follows: 

/Cl,  1-2  = FUNCTION  FCAP(TIN,  TC3,  TOUT,  QCV) 

or 

fu , i:'5~mrrrrTR7'TC3;  Tomv'Qn  j 
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The  function  subroutine  shown  below  must  then  be  entered  as  part 
of  the  input  deck.  The  variables  A,  B,  C and  D appearing  in  the  paren- 
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thesis  correspond  to  TIN,  TC3,  TOUT  and  QC1 , respectively. 

FUNCTION  FCAP  (A,  B,  C,  D) 

IF  (D.GT.l.E-3)  GO  TO  10 
FCAP  = (A-BJ/1.E-6 
RETURN 

10  FCAP  = ( C-B )/l E-6 

RETURN 

END 

2.1.5  Defined  Parameters  Input  Format 

Every  DEFINED  PARAMETER  utilized  in  the  element  input  section  must  be 
defined.  The  parameter  definition  entries  must  be  proceeded  by  the  sub- 
heading card 

/DEFINED  PARAMETERS 

The  format  utilized  to  define  each  defined  parameter  used  is  simply 

/PX  = VALUE 

where  X is  up  to  five  alphanumeric  characters  naming  the  parameter  and 
VALUE  can  have  the  same  formats  used  to  define  elements  except  for  special 
values.  Defined  parameter  values  can  therefore  be  written  in  the  form  of 

• Constants 

• Defined  Parameters 

• Tables 

• Equations 

• FORTRAN  Function  Subroutines 

• Expressions 

The  formats  for  each  have  been  explained  In  detail  In  paragraphs  2. 1.4.1 
through  2. 1.4. 7 and  remain  the  same.  It  will,  therefore,  suffice  here  to 
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give  some  samples  of  possible  entries  under  the  DEFINED  PARAMETERS  card  in 
Q the  input  deck. 

/DEFINED  PARAMETERS 

/PI  = 3.E-12  — Constant  Form 

/POUT  = TABLE  A4(TR3)  -Tabular  Form 

/P3  = P4  —Defined  Parameter  Form 

/P4  = Q3(TR2 , QR1 , 10. ) —Equation  Form 

A> 5 = FCAP(TC1,  TC2,  TIME)  -FORTRAN  Function 

Subroutine  Form 

/P6  = EQUATION  5(10.  , 20.,  P4,  P5)  -Equation  Form 
/P7  = XI  (PI  * P2)  —Expression  Form 

2.1.6  Functions 

The  FUNCTIONS  portion  of  the  Input  is  used  to  define  all  TABLES  and 
EQUATIONS  called  for  under  the  ELEMENTS  and  DEFINED  PARAMETERS  subheadings. 
The  input  cards  for  tables  and  equations  are  preceded  by  the  card 

/FUNCTIONS 

Their  input  formats  are  discussed  in  the  next  two  paragraphs. 

2. 1.6.1  Table  Input  Format 

Every  TABLE  referred  to  under  ELEMENTS  or  DEFINED  PARAMETERS  must  be 
explicitly  defined  under  the  FUNCTIONS  subheading.  The  tabular  data  are 
entered  as  coordinate  pairs  separated  by  commas  preceded  by  the  table  name. 
The  independent  variable  must  be  the  first  coordinate  and  the  dependent 
variable  the  second.  If  the  data  require  continuation  onto  another  card, 
the  card  must  end  with  a dependent  variable  value.  The  next  card  then 
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begins  with  the  next  coordinate  pair. 

The  data  pair  points  must  be  supplied  such  that  the  independent 
variable  is  in  increasing  algebraic  order.  No  specific  limit  to  the  number 
of  point  pairs  that  may  constitute  any  table  is  defined.  Only  single  value 
functions  are  allowed.  It  is  permissible  to  supply  two  consecutive  inde- 
pendent-variable values  that  are  equal,  but  which  have  different  dependent 
variable  values.  This  may  be  done  to  produce  step  functions.  The  example 
below  serves  to  Illustrate  various  acceptable  forms  of  table  input. 

Example 

Referring  to  a previous  example,  the  card 

/qT7  1-2  = TABLE  1 (TCI ) 

was  entered  under  ELEMENTS.  The  function  is  redrawn  here  and  several 
valid  inputs  shown  in  Figure  6. 


Figure  6.  Table  T1  Values  as  a Function  of  TCI 
/TABLE  1 
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1,  0,  0,  0,  1,  1,  2,  1,  2,  1. 


, 1.5,  3,  0,  4,  0 


or 

/n 

/-l,  0,  0,  0,  1,  1,  2,  1 
/2,  1.5,  3,  1.5,  3,  0,  4,  0 
or 

/TABLE"  1 

/=T71 0 
/oTo“ 

/T7T 

/57T 

/27TT5 

/3,  1.5 

/3T0 

/4T0" 

When  the  table  values  are  updated  at  each  solution  time  step,  linear 
interpolation  Is  used  between  the  points  supplied.  For  Independent  vari- 
able values  falling  outside  the  range  of  values  supplied  In  a table,  linear 
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extrapolation  is  performed  to  determine  the  correct  table  value,  and 
therefore,  proper  termination  may  be  necessary.  For  example,  in  Figure  6, 
for  any  value  of  TIME  in  excess  of  3 the  dependent  variable  Q1  will  be 
assigned  the  value  zero  because  of  the  use  of  the  point  pair  4.0.  In 
situations  where  the  same  table  serves  to  define  more  than  one  quantity, 
the  user  need  only  explicitly  define  the  table  once.  For  example,  if  the 
two  heat  sources  Q1  and  Q2  were  input  under  ELEMENTS  as 

/Q1 , 1-2  = TABLE  1 (TCI) 

/Q2,  7-8  = TABLE  1 (TC2) 

only  one  TABLE  I need  be  defined  under  FUNCTIONS. 

2. 1.6. 2 Equation  Input  Format 

Each  unique  equation  referred  to  under  ELEMENTS  or  DEFINED  PARAMETERS 
must  be  explicitly  defined  under  the  FUNCTIONS  subheading  in  terms  of  the 
mathematical  expression  it  represents.  The  format  utilized  to  enter  the 
equation  definition  is  as  follows: 

/EQUAT'lbN  X(br  b2»  . . . bn)'=  (MATHEMATICAL  EXPRESSION) 

/QX ( b1 , b2,  . . . , b^T  = (MATHEMATICAL  EGRESSION) 

X is  the  alphanumeric  name  given  to  the  equation  when  referred  to 
under  ELEMENTS  or  DEFINED  PARAMETERS,  b-j  through  b^  represents  a durruiy 
variable  list  each  containing  up  to  six  alphanumeric  characters  and  repre- 
senting the  variables  called  out  In  the  argument  list  under  ELEMENTS  or 
DEFINED  PARAMETERS.  One  restriction  Is  that  the  DUMMY  VARIABLE  NAMES  MUST 
NOT  BEGIN  WITH  THE  LETTERS  I THROUGH  N INCLUSIVE.  The  mathematical  ex- 
pression must  appear  on  the  card  inside  parentheses  and  be  expressed  in 
terms  of  the  dummy  variables  b^  through  bn. 
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The  mathematical  expression  may  be  any  combination  of  allowable 
operations,  functions,  or  variables.  The  following  mathematical  operations 
and  corresponding  symbols  are  allowed. 


OPERATION  SYMBOL 

Exponentiation  ** 

Multiplication  * 

Division  / 

Addition  + 

Subtraction 


The  order  in  which  operations  are  performed  is  Indicated  by  the  order  in 

which  the  operators  are  listed.  The  use  of  parentheses,  to  denote  clearly 

the  intended  mathematical  combination,  is  suggested  to  avoid  ambiguity. 

For  example,  X + Y * Z should  be  written  as  (X+Y)  * Z ifX+  (Y*Z)  is 

not  intended.  All  standard  FORTRAN  functions  may  be  used,  a few  of  the 

most  commonly  used  are  listed  below. 

FUNCTION  SYMBOL 

Square  root  SQRT  (argument) 

Sine  Function  SIN  (argument  In  radians) 

Cosine  Function  COS  (argument  in  radians) 

Exponential  Function  EXP  (argument) 

Arctangent  Function  ATNA  (argument  in  radians) 

Hyperbolic  Tangent  Function  TANH  (argument  in  radians) 

Natural  Logarithm  Function  ALOG  (argument) 

Common  Logarithm  Function  AL0G10  (argument) 

The  argument  of  any  function  may  be  any  allowable  mathematical  expression. 

Mathematical  expressions  are  not  limited  to  72  characters  (one  card)  and 

may  be  continued  on  subsequent  cards,  using  as  many  as  necessary.  In 

addition,  the  user  may  supply  subprogram  functions  that  he  has  written 

himself  (see  Section  2. 1.4. 7). 
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Example 

Referring  to  a previous  Example,  the  resistor  R3  between  nodes  1 
and  2 was  entered  under  ELEMENTS  as: 

/R3,  1-2  = Q1A'(QR3,  TCI , P4) 

where  f(QR3,  TCI,  P4)  = (QR3  - TC1)*P4. 

The  definition  of  equation  Q1A  is  entered  under  the  subheading 
FUNCTIONS  as: 

/Q1A(A,  B,  C)  = ( (A-B)  * C)) 

where  A,  B and  C correspond  to  QR3,  TCI  and  P4  respectively. 

Example 

The  defined  parameter  PI  = TCI  + QR2  - TABLE  2 (TCI)  is  entered  under 
DEFINED  PARAMETERS  as 

/P 1 = 02A(TCT,"  QR2,  TABLE  2(TCl)) 

The  equation  Q2A  is  In  turn  defined  under  FUNCTIONS  as 

/Q2A(A,  B,  C)  = (A  + B - C) 

At  each  solution  step,  the  ordinate  value  of  Table  II  would  replace  the 
dummy  variable  C and  the  computation  (A  + B - C)  carried  out. 

When  more  than  one  element  or  defined  parameter  can  be  defined  by 
the  same  equation,  only  one  equation  need  be  defined  under  FUNCTIONS. 

For  example.  If  the  two  sources  T1  and  T2  are  defined  as: 

T1  = (10*QC1  - TR2)/TIME 

T2  = ( 1 0*QR4  - TC12)/TIME 

The  two  sources  could  have  been  described  under  ELEMENTS  as 

/Tl,  1-2  = QlOo.,  QC1 , TR2,  TIME) 
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/T2,  B-14  = Q1 (10. , QR4,  TC12,  TIME) 


in  which  case  the  definition  of  Q1  need  only  appear  once  under  FUNCTIONS 
as: 


/Q1(X,  Y,  Z,  T)  = ((X  * Y-Z)/T) 


other  equation  definition  inputs  are  shown  below: 


/Q12(XC1 , P12,  THETA)  = (XC1  - P12  * COS  (THETA)) 

/EQUATION  43(EIN,  ZIN)  = (EIN/ZIN) 

/QT1 6 ( A , B,  C)  = (A  **  B - C/12.E-4) 

2.1.7  Initial  Condition  Input 

Unless  initial  conditions  for  state  variables  are  entered  as  input 
data,  all  voltage  and  current  values  are  assumed  to  be  zero  at  the  first 
solution.  Initial  conditions  other  than  zero  often  exist  and  must  be 
Input.  For  example.  In  AC  analysis.  Initial  conditions  are  needed  by 
the  program  to  linearize  all  non-linear  functions  about  their  operating 
points.  Initial  conditions  follow  the  subheading  card 

/INITIAL  CONDITIONS 


and  are  entered  in  the  forms 

/TX  = NUMBER 
/Q)(  = NUMBER 

where  TX  and  QX  represent  the  temperature  and  the  heat  flow  associated 
with  element  X.  Number  represents  the  value  of  the  temperature  or  heat 
flow  and  must  be  a constant  of  the  form  shown  in  Table  II. 

If  all  Initial  conditions  are  zero,  neither  the  subheading  card  nor 
the  data  are  required.  Care  must  be  taken  to  establish  the  proper 
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polarities  for  initial  conditions.  Also,  the  Initial  capacitor  tempera- 
tures are  positive  when  they  are  consistent  with  the  assumed  temperature 
polarity  for  the  capacitor.  Proper  statements  are  Illustrated  in  Figure  7. 

The  assumed  heat  flow  direction  through  capacitor  Cll  of  Figure  7 is 
from  node  4 to  node  5.  If  the  initial  heat  flow  is  in  this  direction,  it 
is  entered  under  the  INITIAL  CONDITIONS  subheading  as  a positive  quantity. 
If,  however,  the  heat  flows  In  the  opposite  direction,  it  Is  preceded  by  a 
negative  sign.  If  the  actual  initial  temperature  polarity  agrees  with  the 
arbitrarily  chosen  reference  direction,  that  initial  temperature  is  entered 
as  positive;  if  not,  it  is  entered  as  negative. 

© 

0 

+ 

dp  cn 

1 

© 

Figure  7.  Temperature  Polarity  and  Heat  Flow  Direction 
2.1.8  Outputs  Requests 

Any  variable  used  In  the  program  may  be  requested  as  an  output 
variable  under  the  subheading 

/OUTPUTS 

The  variables  listed  under  OUTPUTS  will  appear  In  a printed  listing  as  a 
function  of  time.  The  output  variables  may  be  any  of  the  following: 
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• Temperature  or  heat  flow  associated  with  any 
element  such  as  TCI,  QC1 

• Temperature  or  heat  flow  associated  with  any 
source  such  as  T1 , QE1 , Q2,  TJ2 

• Any  element  value  such  as  Cl 7,  R4 

• Any  DEFINED  PARAMETER  such  as  PI  3,  POUT 

• Any  INTERNAL  parameter.  These  are  parameters 
within  the  program  itself. 

A typical  output  request  is  shown  below: 

/OUTPUTS 

/TCI,  QR1,  R3,  TJ4 
/QX14,  PI  6 

The  variables  may  appear  anywhere  on  the  card  between  column  1 and  72 
separated  by  commas.  As  many  variables  as  desired  may  appear  on  one  card 
and  as  many  cards  as  desired  may  be  used.  If  more  than  one  card  is  used, 
however,  the  last  variable  on  any  card  must  be  followed  by  blanks. 

The  output  is  plotted  If  the  following  format  is  used: 

/n-j , n2,  . . nn,  PLOT 

where  n-j  through  nn  are  desired  to  be  plotted  against  time.  The  output 
requests  below  would  result  in 

/TRY,  TCI,  TSUP,  PLOT 

/QC8 

a listing  and  plotting  of  TRY,  TCI  and  TSUP.  The  fourth  quantity,  QC8 
would  be  output  in  printed  form  only.  If  the  word  PLOT  Is  used,  no  other 
dependent  variables  may  follow  It  on  that  card. 

If  a different  independent  variable  is  desired,  the  following  format 
must  be  used: 
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/QC14,  PLOTfTCW 

In  this  case  the  heat  flow  through  capacitor  Cl 4 would  be  printed  as  well 
as  plotted  as  a function  of  the  temperature  across  it. 

Additional  flexibility  is  available  to  permit  the  user  to  attach  a 
different  label  to  any  output  quantity.  Consider  that  elements  Cl  and  R7 
exist  in  a given  network  and  that  the  temperature  across  both  (TCI  and  TR7) 
are  of  interest.  If  the  user  desires  to  rename  them  as  TIN  and  TOUT,  he 
can  request  the  outputs  as  TCI (TIN),  TR7(T0UT),  PLOT.  The  plot,  rename, 
and  choice  of  independent  variable  options  can  be  presented  in  terms  of 
the  general  format  as: 

YQTY(YLABEL),  PLOT ( XQTY ( XLABEL ) ) 

If  no  rename  is  desired,  the  general  format  reduces  to 

/YQTY,  PLOT (XQTY) 

If,  in  addition,  only  time  is  desired  as  the  independent  variable,  this 
can  be  further  reduced  to 

YQTY, PLOT 

If  no  plotted  Information  is  desired,  the  simplest  form  arises  as 

YQTY,  .... 

2.1.9  Run  Control  Input  Data 

The  RUN  CONTROLS  subheading  card  Is  followed  by  Instructions  specify- 
ing the  type  of  analysis  to  be  performed  and  the  general  parameters  to 
control  the  run.  These  parameters  Include  such  Information  as  the  type 
integration  scheme  to  be  used,  time  limits  set  on  the  run,  etc.  For  all 
the  Run  Control  input  data  available,  the  user  Is  referred  to  a manual  of 
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the  SCEPTRE  program.  These  cards  are  not  processed  by  NELSIM  as  they 
apply  strictly  to  the  execution  of  the  problem  on  SCEPTRE  upon  NELSIM 
output.  A few  of  these  commands  are  illustrated  below. 

2. 1 . 9. 1 Stop  Time 

All  transient  runs  must  have  a problem  duration  in  terms  of  time  units 
calculated  or  assumed.  This  is  the  value  indicating  how  far  out  in  time 
the  solutions  are  to  be  calculated.  This  entry  is  of  the  form: 

/STOP  TIME  - NUMBER 

where  NUMBER  is  a constant  in  time  units  chosen. 

2. 1.9. 2 Maximum  Integration  Passes 

Another  type  of  limit  is  available  that  operates  on  the  number  of 
passes  made  by  the  integration  routine.  In  the  absence  of  any  instruction 
f:om  the  user,  the  program  imposes  a limit  of  20,000  passes.  The  user 
may  change  this  limit  by  entering 

/MAX  INTEGRATION  PASSES  = NUMBER 

2. 1.9. 3 Start  Time 

Most  transient  analysis  problems  begin  at  time  zero  and  for  this  case 
no  entry  is  needed.  If,  however,  the  user  desires  to  begin  the  run  at 
other  than  time  = 0,  the  following  card  Is  necessary 

/START  TIME  = NUMBER 

where  NUMBER  Is  a constant  value  In  units  of  time  assigned  to  the  problem. 

2. 1.9. 4 Minimum  Step  Size 

This  quantity  Is  designed  to  terminate  the  run  automatically  whenever 
f time  solution  Increment  is  required  that  Is  smaller  than  the  minimum 
increment.  If  this  quantity  is  not  supplied,  the  program  will  auto- 
matlcally  compute  a minimum  limit  equal  to  1 x 10  times  the  problem 


31 


I 


duration  (STOP  TIME).  If  a specific  minimum  step  size  independent  of  the 
problem  duration  is  desired,  the  format  is 


For  example  if 


/MINIMUM  STEP  SIZE  = NUMBER 


STOP  TIME  = 5 

.9 

is  specified,  the  minimum  size  would  automatically  be  taken  to  be  5 x 10  . 

If  the  user  is  not  satisfied  with  this  minimum  step-size  and  desires  to 
-5 

change  it  to  1 x 10  he  may  do  so  using  the  following  entry 

/MINIMUM  STEP  SIZE  = IE  -5 

2. 1.9. 5 Punched  Output 

Nominally,  the  output  from  NELSIM  is  in  printed  form  and  punched  card 
output  is  suppressed.  To  obtain  SCEPTRE  input  decks  (a  copy  of  the  printed 
output),  the  punch  option  must  be  turned  on  by  inserting  a card  in  the  run 
controls  input  data  as, 

/PUNCHED  OUTPUT 


2.1.10  Mutual  Inductance  Input  Description 

Although  there  is  no  thermal  analogy  for  an  inductor,  Inductors  will 
appear  in  other  simulations  such  as  the  representation  of  an  electro- 
mechanical system.  Additionally,  the  concept  of  mutual  Inductance  can 
appear  in  a simulation  such  as  for  a transformer. 

Mutual  inductance  is  also  entered  under  the  ELEMENTS  subheading  but  has 
a slightly  different  format.  If  coupling  exists  between  inductors  LI  and 
L2,  the  appropriate  entry  must  include  these  elements  in  place  of  the  node 
identification  such  as 

/MX,  LI-12  = 'VALUE 


w 
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In  addition,  a physical  limitation  of  the  principle  of  mutual  inductance 
must  be  observed  in  order  to  have  a physically  realizable  circuit.  That 
is,  since  coefficient  of  coupling,  k,  is  always  less  than  unity  and  by 
definition 


= M 
•J  LI  L2 


< 1 


the  user  should  be  certain  that  M < JTTT27  Stated  in  words,  the  mutual 
inductance  between  any  two  inductors  must  be  less  than  the  square  root 
of  the  product  of  the  self-inductances  of  the  components  between  which 
the  mutual  inductance  exists.  The  sign  of  M is  positive  if  in  a given 
winding  the  Induced  voltage  of  mutual  inductance  acts  in  the  same  direction 
as  the  induced  voltage  of  self-inductance.  If  the  induced  voltage  of  mutual 
inductance  opposes  the  induced  voltage  of  self-inductance  in  a given  wind- 
ing. M is  negative.  The  proper  sign  for  M for  the  assumed  current  dir- 
ections is  illustrated  in  Figure  8. 


Figure  8.  Mutual  Inductance  Polarities 
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2.1.11  NELSIM  Program  Limits 


Table  III  gives  the  program  data  limits  which  have  been  established 
in  the  program.  These  limits  are  considered  adequate  for  most  normal 
circuit  analysis  work.  Approximate  limits  are  also  given  for  the  number 
of  TABLE  and  EQUATION  functions. 

Table  III.  Program  Data  Limits 


Description  of  Data 

Maximum  Number 

Elements 

300 

Nodes 

301 

Defined  Parameters 

100 

Defined  Parameter  Differential  Equations 

100 

Mutual  Inductances 

50 

Arguments  in  Equation  Value  Specification 

50 

Output  Requests 

100 

Supplied  Initial  Conditions 

300 

Equation  Functions  (1  equation  per  card) 

80  (approximately) 

Cards  per  Equation  Function 

20 

Table  Functions 

80  (approximately) 

2.2  MECHANICAL  SYSTEM  INPUT  LANGUAGE 

Entry  of  a mechanical  problem  into  NELSIM  is  very  similar  to  the  entry 

of  an  electrical  network  as  explained  in  Section  2.1.  An  analog  network 

made  up  of  mechanical  elements  connected  between  nodes  must  be  derived  by 

the  user.  These  elements  may  be  defined  as  linear  or  non-linear  quantities 

using  the  six  formats  outlined  in  Section  2.1.  Two  sets  of  mechanical 

elements  are  allowed  as  input  for  mechanical  problems,  one  for  translational 

systems,  the  other  for  rotational  systems.  The  elements,  their  schematic 
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representation,  as  well  as  program  notation,  are  shov/n  in  Table  IV. 


TABLE  IV.  MECHANICAL  STANDARD  ELEMENTS 


Dashpot 


Spri ng 
Stiffness 

Friction 

Force 

Velocity 


Moment  of 
Inertia 


Torque 

Angular 

Velocity 


Program  Notation 
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The  elements  must  be  entered  in  the  format  shown  in  Section  2.1  and  begin 
with  the  letters  indicated  in  Table  IV.  Table  V shows  the  electrical 
analogs  that  are  generated  from  each  mechanical  element  and  the  numerical 
translation. 

TABLE  V.  MECHANICAL  TO  ELECTRICAL  ANALOG  GENERATION 


Translational 

Rotational 

FORCE  *>  CURRENT  (I  = F) 

VELOCITY  =>  VOLTAGE  (V  = V) 

MASS  =>  CAPACITANCE  (C  = M) 

SPRING  =>  INDUCTOR  (L  = 1/K) 
FRICTION  *>  RESISTANCE  (R  = 1/B 

TORQUE  =t»  CURRENT  (I  = T) 

VELOCITY  o VOLTAGE  (V  = W) 

MOMENT  OF  I CAPACITANCE  (C  = J) 

SPRING  =>  INDUCTANCE  (L  = 1/K) 

FRICTION  =>  RESISTANCE  (R  = 1/B) 

2.2.1  Mass  Connections 

The  node  names  are  user  derived  in  general.  For  mass  connections, 
however,  a restriction  Is  placed  on  the  user.  All  masses  must  be  connected 
to  a ground  node  beginning  with  the  letter  G.  Further,  all  masses  must 
appear  at  node  points  only.  If  two  or  more  masses  are  connected  to 
different  points,  each  must  be  connected  to  the  common  ground  node.  The 
following  system  shows  a possible  nodal  connection  for  a 2 mass  system. 
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2.2.2  Friction  Connections 

Friction  can  occur  as  the  resistive  force  between  two  masses  or  be 
represented  as  a dashpot  (Figure  9).  Two  types  of  nodal  connections  are, 
therefore  allowed  for  a friction  input,  one  between  masses  and  one  between 
nodes. 


© © 


Figure  9.  Frictional  Elements 
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The  frictional  element  inputs  for  Figure  9 are  shown  below: 


ELEMENTS 


B1 , Ml -M2  = VALUE 
B2,  2-G  = VALUE 
B3,  M2-G  = VALUE 

2.3  MECHANICAL  SYSTEM  STATE  VARIABLES 

The  state  variables  applicable  to  mechanical  system  and  recognized  by 
the  program  are  shown  in  Table  V.  These  can  be  referred  to  in  the  input 
language  by  placing  the  designation  letter  of  Table  V before  the  element 
it  applies  to,  for  example  FK1  represents  the  force  exerted  on  spring  K1 . 
TABLE  VI.  MECHANICAL  SYSTEM  STATE  VARIABLES 


Translational 

Rotational 

State 

Variable 

Program 

Notation 

State 

Variable 

Program 

Notation 

Force 

F 

Torque 

T 

Velocity 

V 

Angular 

Velocity 

W 

The  state  variables  are  analogous  to  voltage  and  current  in  the  electrical 
case  where  angular  velocity  represents  voltage  and  torque  represents 
current. 

INITIAL  CONDITIONS 

Initial  conditions  are  handled  In  a manner  similar  to  the  electrical 
case.  The  mechanical  analog  of  a capacitor  Is  a mass  (M)  for  the  transla- 
tional case  and  the  moment  of  Inertia  (J)  for  the  translational  case.  The 
mechanical  analog  of  an  Inductor  Is  the  spring  stiffness  K.  Therefore,  the 
Initial  conditions  of  a mechanical  system  can  be  entered  as  the  force 
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exerted  on  Ks  and  velocity  on  Ms  and  about  Js.  Some  samples  are  given 
below. 

INITIAL  CONDITIONS 

WJ5  = 10 

WM1 0 = 15 

VM2  = 5 

FK1 2 = -20 

2.4  OUTPUTS 

The  outputs  are  entered  in  the  manner  described  in  Section  2.1.10  with 
the  electrical  element  names  replaced  by  mechanical  element  names,  voltage 
replaced  by  velocity  (V)  or  angular  velocity  (W)  and  current  replaced  by 
force  (F)  or  torque  (T). 

Mechanical  System  Input  Examples.  Consider  the  mechanical  problem 
shown  in  Figure  10. 


Figure  10.  Mechanical  Problem 

The  noded  mechanical  system  and  input  for  NELSIM  is  shown  in  Figure  11. 
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ANALOG  DESCRIPTION 
MECHANICAL 
ELEMENTS 
VI ,6-3=Ql (TIME) 

K1 ,62-3=1 
K2,61 -2=P2 
K3, 6-1=3 
M2, 2-62=. 2 
M3, 1-61=. 3 
B1 ,61-2=1 
FUNCTIONS 

Q1(T)=  ( COS( 2 . *T) ) 

DEFINED  PARAMETERS 
P2=2 
OUTPUTS 

VM2(V2),VM3(V3)  .PLOT 
RUN  CONTROLS 

INTEGRATION  ROUTI NE= I MPL I Cl T 
STOP  TIME =10 

END 


Figure  11.  NELSIM  Input  for  Mechanical  System 


Figure  12  represents  a rotational  system  and  its  input  for  the 


normalized  case  where  all  values  are  set  equal  to  1. 


ANALOG  DESCRIPTION 

MECHANICAL 

ELEMENTS 

Wl,  G-l  = 

1 

Kl,  1-2  = 

1 

J2,  2-G1  = 

1 

K2,  Gl-3  = 

1 

J3 , 3-G2  = 

1 

B3,  J3-G  = 

1 

B2,  J2-G  = 

1 

OUTPUTS 

WJ2 , WB2 , 

PL0T 

RUN  CONTROLS 

STOP  TIME 

= 1 

END 

K1  «-2(t)  *2  u’3(t) 


Figure  12.  NELSIM  Input  for  a Rotational  Problem 
2.5  THERMAL  SYSTEM  INPUT 

The  Input  to  NELSIM  for  a thermal  problem  is  similar  to  the  input  of 
a mechanical  problem.  The  elements  for  a thermal  system  are  shown  In 
Table  VII. 
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TABLE  VII.  THERMAL  ELEMENTS  AND  STATE  VARIABLES 


Element 

Schematic  Representation 

Program  Notation 

Thermal 
Capaci tance 

1( 

C 

Thermal 

Resistance 

Wv 

R 

Temperature 

Source 

T 

Heat  Flow 
Source 

Q 

iSBi 

The  user  must  derive  a thermal  analog  network  made  up  of  the  elements  shown 
In  Table  VII.  The  value  of  these  elements  may  be  entered  in  any  of  the 
formats  documented  in  Section  2.1.6. 

Thermal  State  Variables 

The  state  variables  for  thermal  input  are  Temperature  (T)  and  heat 
flow  (Q)  and  are  analogous  to  voltage  (V)  and  current  (I)  respectively  in  the 
electrical  network. 

Initial  Conditions 

Initial  conditions  for  a thermal  system  may  be  entered  as  the  temper- 
ature across  a thermal  capacitance.  This  is  analogous  to  a voltage  across 
an  electrical  capacitance.  Some  examples  follow. 
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INITIAL  CONDITIONS 
TCI  = 10 
TC2  = 20 

Outputs 

Outputs  are  requested  as  defined  in  Section  2.1.10  with  the  electrical 
elements  replaced  by  the  thermal  elements  of  Table  VII,  voltage  (V) 
replaced  by  Temperature  (T)  and  current  replaced  by  heat  flow  (Q). 

Thermal  Input  Example 

The  following  problem  Illustrates  the  Input  of  a thermal  system. 
Consider  the  rocket  engine  test  setup  shown  in  Figure  13  for  which  it  is 
desired  to  study  the  transient  temperature  response  within  the  engine  wall. 
The  input  of  the  problem  requires  the  various  subdivisions  within  the  wall 
to  be  indicated  by  nodes  (1-3)  and  shown  in  Figure  13  (b).  The  thermal 
resistance  between  subdivisions  is  simulated  using  resistors  and  thermal 
capacitance  using  capacitors.  The  element  input  is  shown  below. 


ANALOG  DECRiPT  ION 
THERMAL 
ELEMENTS 
TIN,G-l=Tl 
Kl, 1-2=1 
C1,2-G=J 
R2, 2-3=1 
C2,3-G=l 
R3, 3-4=1 
TUUT,G-A=1 
FUNCT I DNS 

T l=Ot 1 , 1,1,  2,10,  S,1U*  6,1,  10,1 

OUTPUTS 

TIN,QC1,QC2,QR3.PL0T 
TR3, TCI, PLOT 
RUN  CONTROLS 

INTEGRATION  ROUT  1 NE  = 1 MPL 1 C 1 T 
STOP  TIME=1S 
RUN  INITIAL  CONDITIONS 
END 
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The  elements  have  all  been  entered  as  constant  values  for  simplicity  but 
in  general  will  be  non-linear  functions  of  the  materials  involved  and  will 

be  entered  using  a variety  of  value  formats. 

Test  Stand 


Plane  P 


Figure  13.  Heat  Transfer  Example 
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2.6  ELECTRO-MECHANICAL  ELEMENTS  AND  MODELS 

The  entry  of  electro-mechanical  problems  into  NELSIM  is  possible  in 
several  ways.  The  differential  equation  and  transfer  function  entries  are 
defered  until  a general  treatment  of  these  topics  in  Sections  III  and  IV, 
respectively.  Two  additional  modes  of  entry,  however,  are  available  to  the 
user,  namely 

1)  Electro-mechanical  element  entry 

2)  Model  entry 

2.6.1  Electro-Mechanical  Element  Entry 

NELSIM  will  allow  for  the  input  of  an  electro/mechanical  system  in 
terms  of  the  lumped  elements  of  both  the  electrical  and  mechanical  portions 
of  the  system  coupled  together  through  common  defined  parameter  names, 
algebraic  or  first  order  differential  equations  or  dependent  sources.  The 
mechanical  and  electrical  elements  are  entered  folio,  ng  the  subheadings 
ELEMENTS  (MECHANICAL)  and  ELEMENTS  (ELECTRICAL),  respectively.  The 
electrical  elements  allowed  are  discussed  in  Section  2.1.  The  input  is 
better  illustrated  through  an  example.  Consider  the  galvenometer  sketch 
provided  In  Figure  14.  A galvanometer  ideally  produces  a deflection  that 
is  proportional  to  an  electrical  input. 


Figure  14.  Galvanometer  Schematic  Representation 
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Figure  15  represents  the  equivalent  lumped  network  in  terms  of  its 
electrical  and  mechanical  parts 


Figure  15.  Lumped  Element  Galvanometer  Equivalent 

where  J and  B represent  the  moment  of  inertia  of  the  rotor  and  the  viscous 
damping  that  results  from  air  friction.  The  element  KR  represents  the 
rotational  stiffness  produced  by  a spring  attached  to  the  rotor.  The 
constant  D is  defined  to  be  fit  a where  £is  the  magnetic  flux  density,  ( is 
the  total  length  of  the  conductor  and  a Is  the  radius  of  the  rotor. 

The  Input  language  where  all  the  elements  are  normalized  to  one 
follows. 


ANALOG  DESCRIPTION 
ELECTRO  MECHANICAL 
ELEMENTS  ( ELECTRICAL'' 


EltG-1-1  Electrical  Source  in  terms 

Rl»l-£-l  of  rotational  Velocity 

EC  f6-2*QEC  <PD  ! ‘ 


ELEMENTS < MECHANICAL > 

Tl  •G-.3*=OEC»  PD » IRl  Applied  Mechanical 

JR»3-G*1  Torque  in  terms  of 

BR>3-G*1  current  through 

KR  *3-G*  1 resistor 
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The  two  dependent  sources  EC  and  T1  provide  the  interface  between  the  two 
systems  and  are  entered  directly  in  terms  of  the  elements  involved. 

Element  NAME  uniqueness 

NELSIM  converts  all  mechanical  elements  to  electrical  equivalent 
elements  by  replacing  the  first  letter  of  the  element.  For  example  Mass 
M15  becomes  capacitor  C15  and  spring  constant  K1  becomes  inductor  LI.  Care 
must,  therefore,  be  taken  by  the  user  to  avoid  duplication  of  names  in  the 
resulting  network.  This  can  be  accomplished  by  checking  that  all  electrical 
and  mechanical  elements  have  unique  denotations  following  the  first  letter. 
2.6.2  Built  in  Models 

This  option  will  allow  the  user  to  include  models  of  generally  used 
electro-mechanical  devices  in  his  system  by  simply  calling  for  them.  It 
is  difficult  to  derive  models  general  enough  to  suit  all  needs  but  it 
should  be  kept  In  mind  that  the  capability  to  input  system  problems  in 
terms  of  differential  equations  or  transfer  function  blocks  allows  the 
user  to  build  his  own  models.  The  two  models  presently  implemented  into 
NELSIM  are  that  of  a rate  gyroscope  and  an  accelerometer.  The  call  to  a 
model  will  be  of  the  form 

T1 , N1-N2  = MODEL  XXXX 

where  N1  represents  the  input  to  the  model  and  N2  the  output  and  XXX  is 
the  built  in  model  name.  The  two  models  are  discussed  below.  The  model 
entries  are  submitted  following  the  subheading  MODELS. 

2 . f . 2 . 1 Rate  Gyro  Model 

Kate  gyros  are  used  widely  in  any  application  requiring  tracking  and 
stabilization.  Aircraft  and  missiles  utilizes  gyros  for  stabilization 
when  angular  deviations  are  noticed.  Rate  gyros  are  also  used  extensively 
in  conjunction  with  tracking  antennas  so  that  the  angular  velocity  of  the 
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tracked  vehicle  can  be  measured.  A mathematical  block  diagram  of  a rate 
gyro  is  given  in  Figure  16. 


Figure  16.  Block  Diagram  of  Rate  Gyro 

In  the  figure  w..  represents  the  rate  input  to  the  gyro  and  eQ  the  pro- 
portional output  voltage  generated.  The  K's  represent  amplifiers  within 
the  system  to  attain  the  desired  signal  levels  with  KSG  representing  the 
signal  generator  which  generates  an  electrical  signal  proportional  to 
the  displacement  angle  0Q.  The  transfer  function  between  12  and  T1  is 
calculated  as  follows: 


T2 

T1 


S(Hf  S) 

iTlZc- 
so+j-  S) 


when  simplified  the  above  equation  becomes 


T2 

TT 


1 + 


2 
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The  above  represents  the  closed  loop  transfer  function  of  the  rate  gyro 
which  can  be  recognized  as  having  the  form 


F(S)  = 


+ 2i 


S + 


where  4 is  the  damping  factor  and  u;n  the  natural  frequency. 
Like  terms  yield  that 


j-  and  £ = 


U KJ 


which  allows  the  user  to  adjust  both  of  these  quantities  by  adjusting 
' the  various  gains. 

The  model  equivalent  for  NELSIM  input  is  shown  in  Figure  17  where 
PB  = 1/C  and  PC  = J.  The  feedback  gain  PF  has  been  added  for  further 
generalization  of  the  model  and  if  set  equal  to  one  will  result  in  the 
transfer  function  F(s).  Use  of  the  model  requires  the  user  to  input  the 
request  under  the  ELECTRO-MECHANICAL  heading  and  following  the  MODELS  sub- 
heading card.  Usage  also  requires  the  user  to  define  all  the  gains  needed 
by  the  model  (PA,  PB,  PC,  PD,  PF).  The  letter  P is  attached  to  both  the 
input  node  and  the  output  node  of  a model  as  well  as  the  name  assigned  by 
the  user.  If  the  model  is  referred  to  as  Ml  for  example  the  input  to  the 
model  becomes  PINM1  and  output  P0UTM1 . The  model  name  Is  also  attached  to 
the  internal  parameters  to  maintain  parameters  name  uniqueness  ( 1 . e . , 

PIM1 , PG2M1 , PFM1 ) . To  illustrate  consider  Figure  18  for  which  H 
represents  the  gyro  model . 
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Figure  17.  NELSIM  Equivalent  Rate  Gyro 


Figure  18.  NELSIM  Gyro  Model 
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The  input  to  the  NELSIM  program  is  shown  below. 


ANALUb  UttklFTlUN 
LLECTkO  litCHAMCAL 
FiUUtLS 


1 I , i-2=uYkti 


Input  Step — — 

Step  

Description 


ULt  iKLU  FARAMLTEKS 
FAT  1-1 
Fbl 1-1 


FU1  1 = 
FC1  1 = 
FF  i 1- 


Model  Name  Added 


to  Parameters 


F1T1-T1 
FUkCT  iUUS 

T1=0,0,  1,0,  2,1,  10,1 

UOTFUTS 

FIT  1 ( INFUT).F2T 1 (OUT ) , Fu2T 1 ,FLUT 
k(jh  CUkiKULS 


Response  request 
output  node 


at 


i>  I LiF  i 1 Fit  — l U 

IkTLUkAT ION  kOUT INE-lMFLiOll 
tNU 


2 . 6 . 2 . 2 Accelerometer  Model 

The  accelerometer  is  the  basic  measuring  element  In  an  inertial  navi- 
gation and  guidance  system.  It  is  the  function  of  the  accelerometer  to 
provide  measurement  and  Integration  of  linear  acceleration  to  successfully 
accomplish  the  functions  of  position  and  velocity  Indication  and  generation 
of  proper  steering  signals.  Although  numerous  models  are  available  the 
double  integrating  accelerometer  shown  in  Figure  19  has  been  implemented 
into  NELSIM. 

In  the  model  a^  represents  the  input  acceleration,  eQ  the  signal 
voltage  generated  and  9 the  angular  acceleration  of  the  motor  shaft 
resulting  from  the  application  of  eQ.  The  double  Integration  of  the 
angular  velocity  (9)  yields  the  angular  velocity^  and  relative  angle 
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(9)  of  the  shaft.  These  two  quantities  are  proportional  to  the  craft 
velocity  (V)  and  distance  traveled  (X).  The  mathematical  formulation  is 
given  below. 


3 


where: 


T3  _ ^K2KSGK4^S 

TZ  S(1  + ^3  S)  (1  + K5S)  + (K2K$gK4)S 
K2 

= N 

1 + aS  + bS2  + N 


N = (K2KsgK4)S 


inserting  T = and  applying  the  final  value  theorem  for  an  assured 
step  input  yields 


V 


KVi 


in  the  steady  state  which  shows  the  angular  acceleration  of  the  shaft  to 
be  proportional  to  the  Input  acceleration.  Another  relationship  exist 
for  T3,  namely 

T3  = Kg9S2 

equating  the  two  yields 

KKlai  = K60s2 

To  determine  velocity  and  distance  changes  as  a function  of  a.  both  sides 
of  this  equation  are  integrated. 


S3 


A9dt 


The  model  programmed  Into  NELSIM  is  illustrated  in  Figure  20. 


Figure  20.  NELSIM  Accelerometer  Equivalent 
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The  usage  and  response  of  the  model  is  illustrated  for  the  case  where  all 
gains  have  been  normalized  to  1 and  the  input  to  the  model  is  a unit  step. 
Consider  Figure  21  where  H is  the  accelerometer 


IN 

O 


Figure  21.  NELSIM  Accelerometer  Model 


The  input  to  NELSIM  is  given  below: 


ANLAEOO  UtCklPT ION 
LLEUkU  HLCHAWiCAL 
MODELS 

l 1 , iN-UU=ACCELEROME  I Ek 

OEF  ii'tLU  PARAMETERS 

HAT  1=1 

PbT  1 = 1 

PCT  1 = 1 

POT  1 = 1 

PEI  1 = 1 

PFT 1=1 

POT  1 = 1 

P I NT  1 =T 1 

F UNC  I I DNS 

T1=U,U,  0,1,  10,1 

OUTPUTS 

P I NT  1 , PAT  1 , PA2 T 1 , PA3  T 1 , PAA T 1 , PLOT 
PAS  1 1 ,PUU1 1 .PLOT 
KUh  CONTKULS 

INlEbRAT  ION  KUUT  INE=IMPL ICi  1 

Slop  T1HE=10 

END 


As  was  required  for  the  gyro  model,  all  gains  within  the  model  must  be 
provided  by  the  user  with  the  model  label  (T1 ) added  to  each  name. 

2.7  ELECTRO-OPTICAL  SYSTEM  INPUT 

Electro-optical  systems  are  represented  in  NELSIM  as  equivalent 
electrical  circuits  capable  of  simulating  the  AC  characteristics  of  photo 
diode  devices. 


J 
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The  most  commonly  used  circuit  is  illustrated  in  Figure  22  followed 
by  the  NELSIM  input  describing  the  circuitry. 


Figure  22.  Equivalent  Circuit  of  a Photo  Diode 


ANALOG  L/LSCixIPl  I UN 
ELECTRO  OPTICAL 
ELEMENTS (ELECTRICAL) 

JPN,G-  l=OEO(PEFi  ,PCHG,PCON,PUER,  TPL) 

K 1 , 2-G-6 

Cl, 3- G=3 

US, 3-4=12 

f<L,  4-0=12 

FUliCT  IONS 

TPL=0,C,.2,2,.4,6,.05, 10,. 6, 12,0, 13 
GEO(A,B,C,L>,E)  = (A*B/(C*D)*E) 

REPINED  PARAMETERS 

PLFF=. 25 

PCI1G= 1 E- 7 

PC0N=1 23L-4 

POf  R— 3 7 

UOTPUTS 

JPU, JKL,PLUT 

4 . 0 1 . CUN  l i\bLS 

INI  LGKA I 1 UN  RUG  l I NL—  IMP  1C  1 l 

STUP  1 1/  iL=  1 3 

EM) 
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In  the  above  input,  the  photo  current  I PH  is  described  by  a perfect-square 
law  formulation 


IpH  = (n  q P/hv) 

where, 

’l  = quantum  efficiency 

P/hv  = average  number  of  incident  photons 
per  unit  time 

IpH/q  = average  number  fps  unit  time  of 

electrons  emitted  from  the  photocathode. 

Since  the  input  is  as  an  equivalent  circuit,  the  circuit  elements 
are  detailed  as  ELEMENTS  (ELECTRICAL)  and  the  element  names  are  consistent 
with  the  network  preparation  rules  of  Section  2.1.2. 

Since  the  device  response  is  dependent  upon  the  frequency  of  the 
energy  source  it  is  best  to  construct  a circuit  for  frequency  bands  and 
sum  the  produced  currents  linearly,  or  to  restrict  the  analysis  to  a 
certain  frequency  band. 

2.8  UNITS  AND  SCALING 

The  following  paragraphs  describe  the  scaling  option  in  NELSIM.  In 
order  to  allow  NELSIM  to  be  completely  compatible  with  SCEPTRE  In  the 
flexibility  of  its  input  language  the  scaling  module  Is  semiautomated  and 
requires  user  interaction. 

2.8.1  Units 

SCEPTRE  allows  any  set  of  parameter  units  as  long  as  the  element 
magnitudes  with  respect  to  these  are  consistent.  To  maintain  the  two  input 
languages  consistent,  the  same  will  be  true  of  NELSIM.  The  user  is  allowed 
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to  enter  a problem  in  any  set  of  units  he  chooses  and  is  responsible  for 
compatibility  throughout  his  system. 

2.8.2  Scaling  Guidelines  and  Automation 

One  of  the  problems  encountered  in  system  analysis  utilizing  numeri- 
cal integration  is  that  of  maintaining  a consistent  set  of  parameter  units 
such  that  the  magnitudes  of  all  the  solution  variables  are  consistent 
with  the  specified  error  criteria.  Numerical  Inaccuracies  result  and 
excessive  solution  time  is  spent  when  units  are  not  carefully  chosen  due 
to  the  difficulties  in  meeting  specified  error  criteria.  Previous 
experience  in  the  system  analysis  area  and  application  of  SCEPTRE  and 
other  codes  to  large-scale  circuit  and  system  problems  has  shown  that  is 

desirable  to  maintain  the  magnitude  of  the  solution  variables  in  the  range 
-3  3 

between  10  and  10  . For  example,  voltage  and  current  are  the  two 
solution  variables  used  by  SCEPTRE  and  It  is  desirable  to  keep  these 
variables  within  the  limits  shown  below 

10'3<  I V I < 103 
l(f3S  I I I HO3 

The  desired  magnitude  ranges  can  be  accomplished  through  knowledge  of  the 
basic  laws  governing  a particular  system  and  a knowledge  of  the  ranges 
expected.  If  the  expected  units  of  the  solution  variables  and  time  are 
known,  the  system  element  units  can  be  scaled  to  yield  those  units.  The 
above  is  exemplified  through  a high  speed  transistorized  circuit  in  which 
the  response  is  known  to  be  in  terms  of  volts,  milliamps  and  nanoseconds. 
Equations  1 through  3 Illustrates  the  applicable  relationships  between 
state  variables  and  electrical  elements. 


58 


(1) 


I 


c ■ I f (2) 

l ■ v ji  (3) 

Where:  V = Voltage 

I = Current 
R = Resistance 
C = Capacitance 
L = Inductance 


To  obtain  the  response  in  terms  of  volts,  mllliamps  and  nanoseconds  the 
elements  are  scaled  as  follows: 


R = T = = = ]°3  (kilohms) 

I ma  10"^ 

C = I = ma,  5§ff-s  = ^1n"  ) = 10"12(picofarads) 

L = V = Volts.  ^ ^ = 10"6(microhenries) 

□ I md  1 q-j 

Definition  of  the  network  elements  in  terms  of  these  units  will  bring  the 
solution  variables  closer  to  their  desired  range. 

The  same  procedure  is  applicable  to  other  systems.  Equations  4 through 
8 present  the  relationships  necessary  for  mechanical  and  thermal  systems. 
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Where: 


MECHANICAL  SYSTEMS 


(4) 


M = F % (5) 


, _ 1 dF 
K “ V dt 


(6) 


V = Velocity 
F = Force 
M = Mass 

K = Spring  Constant 
B = Viscocity/Friction 

THERMAL  SYSTEM 

c ■ i af  <7> 


R = 


(8) 


I 


Where:  Q = Rate  of  heat  flow 

C = Thermal  Capacitance 
R = Thermal  Resistance 
T = Temperature 

2.8.3  Automated  Scaling 

An  algorithm  has  been  implemented  within  the  NELSIM  program  to  achieve 
scaling  of  a network.  The  program  generates  an  electrical  analog  from 
various  types  of  system  inout.s.  The  scaling  algorithm  scales  all  the  ele- 
ments to  magnitudes  consistent  with  user  provided  guidelines.  The  option 
require  the  user  to  input  the  expected  units  of  the  solution  variable  and 
time  in  terms  of  electrical  quantities.  Table  VIII  shows  the  analogies 
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between  the  various  systems  involved.  Both  mechanical  and  thermal  quanti- 
ties are  translated  into  the  electrical  quantities  shown  in  the  table, 
such  that  in  the  translated  network  the  only  solution  variables  are 
electrical.  The  user  will  enter  a problem  in  terms  of  the  standard  units 
familiar  to  the  field  involved.  Table  IX  illustrates  the  units  that 
are  expected  if  the  scale  option  is  utilized.  The  user  indicates  the 
scaling  to  be  performed  by  entering  the  expected  magnitudes  of  the 
solution  variables  and  time.  For  example,  consider  the  electronic  network 
discussed  earlier  entered  in  terms  of  the  standard  elements  shown  in 
Table  IX.  To  obtain  the  solution  in  terms  of  volts,  milliamps  and  nano- 
seconds the  network  must  be  scaled.  The  user  input  would  consist  of  the 
card  shown  below  which  states  that  the  scaling  option  is  desired 

SCALE  OPTION,  V0LTAGE=1 , CURRENT=10E-3,  TIME=10E-9 

and  the  voltage  is  to  remain  in  terms  of  volts  but  current  and  time  are  to 
be  scaled  to  milliamps  and  nanoseconds  accordingly.  Use  of  the  scaling 
option  for  all  system  components  will  maintain  a consistent  set  of 
variable  values. 

The  output  of  the  program  consists  of  two  listings.  The  first  is  a 
listing  of  the  analog  network  generated  with  no  scaling  applied.  The 
second  consists  of  the  analog  network  with  all  constant  elements  scaled 
and  a list  of  the  scale  factors  used  on  the  elements.  The  user  is  then 
required  to  scale  all  non-constant  element  values  to  complete  the  scaling 
task. 

To  illustrate  the  use  of  the  scaling  option  consider  the  mechanical 

prob' em  presented  earlier.  It  is  desired  to  scale  time  by  a factor  of 
-3  -3 

10  , force  by  a factor  of  10  and  leave  the  velocity  units  the  same. 


61 


TABLE  VIII.  SYSTEM  ANALOGIES 


I 


SYSTEM 

FUNCTION 

ELECTRICAL 

MECHANICAL 

THERMAL 

Translational 

Rotational 

SOLUTION 

VARIABLES 

Voltage 

Velocity 

Angular  Velocity 

Temperature 

Current 

Force 

Torque 

Heat  Mo 

TIME 

Time 

Time 

T Ime 

T Ime 

ELEMENTS 

Resistance 

Friction 

Friction 

Thermal 

Resistance 

Capacitance 

Mass 

Moment  of  Inertia 

Therral 

Capacitance 

Inductance 

Spring  Stiffness 

Spring  Stiffness 

I TABLE  IX.  STANDARD  SYSTEM  UNITS 


UNITS 

PHYSICAL  MEDIUM 

ELEMENTS 

SYMBOL 

STANDARD 

MKS 

ENCLISH 

rv 

Electrical 

Resistance 

m 

Ohms 

Capacitance 

Farads 

Inductance 

Hi 

Henries 

Mechanical 

Mass 

M 

Kilogram 

Slug 

ii 

Spring 

Stiffness 

K 

Newton/ 

meter 

lb/ ft 

Friction 

B 

n-sec/m 

lb/ft 

Moment  of 
Inertia 

j 

ml  . 

(rad/sec*) 

ft  lb/  - 
( rad/sec*) 

■ 

Angular 

Velocity 

W 

rad/ sec 

raCs/sec 

■ 
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TABLE  IX.  STANDARD  SYSTEM  UNITS  (CONTINUED) 


PHYSICAL  MEDIUM 

ELEMENTS 

SYMBOL 

STANDARD 

MKS 

ENGLISH 

BTU 

Thermal 

Temperature 

T 

°F 

Heat  flow 

rate 

Q 

BTU/sec 

Capacitance 

c 

BTU/°F 

Resistance 

R 

(BTU/°F/ 

lb 
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As  shown  in  Table  VIII,  the  analog  quantities  for  force  and  velocity  are 
current  and  voltage,  respectively.  The  scale  request  is 

SCALE  OPTION,  TIME=lE-3,  CURRENTS E-3 , V0LTAGE=1 

The  output  of  NELSIM  will  first  list  the  network  analog  generated 
with  no  scaling  applied,  and  the  second  lists,  the  value  of  analog 
elements  with  all  constant  elements. 
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i SECTION  III 

TRANSFER  FUNCTION  INPUT  LANGUAGE 

The  input  language  for  a network  topology  was  described  in  the  pre- 
ceding section.  For  a certain  class  of  problems,  a signal  flowgraph 
representation  of  the  system  may  already  be  at  hand,  for  example,  control 
systems . 

An  input  language  has  been  devised  to  allow  the  direct  input  of  a 
flowgraph  into  NELSIM.  The  program  interprets  this  input,  and  from  it 
derives  the  set  of  first  order  differential  equations  necessary  to  execute 
the  problem  in  SCEPTRE. 

The  input  utilizes  the  same  type  of  format  available  for  circuit 
inputs  in  SCEPTRE  and  is  nearly  identical  to  SYNAP.  The  data  are  placed 
under  headings  and  subheadings  in  a similar  fashion  and  is  format  free. 

The  statements  consist  of  descriptive  user-derived  names  for  nodes,  com- 
ponents, and  value  specifications  separated  by  delimiters,  such  as  contnas , 
parentheses,  and  dashes.  The  input  data  may  be  punched  anywhere  on  the 
card  within  the  first  72  columns.  To  continue  a statement  from  one  card 
to  another,  the  only  requirement  is  that  the  discontinuation  be  made 
irrmedi ately  after  delimiters,  with  the  delimiter  being  the  last  non-blank 
character  on  the  card.  The  input  data  organization  is  shown  in  Figure  23. 
Each  entry  is  described  in  the  following  paragraphs. 
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Main 

Headings 


- TRANSFER  FUNCTIONS 


Comments 

Brandies 


DEFINED  PARAMETERS 


Run  Controls 
Outputs 
Sensi  ti  vities 


> Subheadings 


— RERUN  DESCRIPTION  (N) 


Values 

Run  Controls 


— END 


| Subheadings 


Figure  23.  Graph  Description  Data  Organization 
3.1  FLOWGRAPH  INPUT  PREPARATION 


The  flowgraph  description  is  defined  in  terms  of  branch  transmi ttances 
connected  between  nodes.  The  steps  listed  below  should  be  followed  in  pre- 


pari ng 


l 

■ 


a signal  flowgraph  for  input. 

Assign  a name  or  number  to  all  nodes  in  the  graph 
(four  or  less  alphanumeric  characters). 


Assign  a name  or  value  to  all  transmittance  parameters 
(four  or  less  alphanumeric  characters).  All  parameters 
must  begin  with  the  letter  P and  be  defined  under  the 
DEFINED  PARAMETERS  subheading. 


Label  signal  flow  direction  between  node  connections 
using  arrowheads 

Figure  24  illustrates  a flowgraph  ready  for  input  generation. 
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3.2  GRAPH  DESCRIPTION 


I 


For  Network  topologies  the  first  card  was  ANALOG  DESCRIPTION,  in  the 
same  matter,  the  first  card  to  describe  a signal  flowgraph  is 

/transfer  FUNCTIONS 

This  card  indicates  to  the  program  that  a flowgraph  is  being  input  and 
allows  the  program  to  enter  the  graph  processor  directly. 

3.3  BRANCHES 

The  flowgraph  is  made  up  of  branches  containing  some  transmittance 
characteristics  between  one  node  and  another.  The  format  for  any  branch  is 
shown  below  where: 

/NODE 1 - N0DE2  (BRANCH  TRANSMITTANCE) 

• N0DE1  denotes  the  initial  node  from  which  the  signal 
originates  (from  node) 

• N0DE2  denotes  the  final  node  of  the  transmittance  ( to 
node) 

The  general  form  of  the  branch  transmittance  is  a rational  polynomial  of 
the  form 


p(s)  = 

N0*B01*B02  * • * 

B0n  + N1  *B11  *B1 2 ‘ 

••BlnS+ 

...N  *B  , *B  , *B  „...B 
m ml  ml  m2  i 

M0*C01  *C02' * 

•con + Mrcirci2-- 

•Cln  S + ■ 

k kl  k2  kn 

where 

Nq,  N] , N2 
B01 ’ B02’ 

...Nm,  M0,  MP  . 

‘ • ,B0n  ’ B1 1 * B1 2 ’ • 

" *Mk 
’ ,Bln  ’ 

are  transmittance 
coefficient  parameters 

C c 

01  ’ ^02  ’ 

...c0n,  cn,  c12  . 

• -Cln 

Polynomials  of  this  form  are  recognized  by  the  program  by  using  the  following 
input  form: 
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/NODE  1 - NODE2  (NUMERATOR/DENOMINATOR) 
where  the  expanded  form  yields 

/NUl)t  - NODE  (nq*B01*IJ02  BOn’  C0’  K1  *B]  1 *^1 2 ~ Bln>  ^ * 


/N7B^2  ...  B2n,  2 ...  Bmn,  »/ 


/M~*r  *r 

1 o 01  02 


^On  ’ •••  ^in>  M9*C91  *C99  ... 


1 11  12 


In  1 


'2  21  22  •••  °2n’ 


I 


/...,  Mk-Bkl*Bk2  ...  Bk^KT 


3.4  INPUT/OUTPUT  REFERENCE 

Each  node  name  entered  in  the  branch  description  is  translated  into  a 
defined  parameter  by  NELSIM  by  adding  the  letter  P to  each  node.  This  is 
done  to  satisfy  the  SCEPTRE  input  as  well  as  interfaces  requirements  that 
will  be  discussed  later.  The  user  desiring  to  define  or  output  a node  must 
refer  to  it  as  PNODE  in  the  input.  An  example  will  illustrate  the  procedure. 
The  branches  of  the  feedback  system  shown  in  Figure  25  are  described  as 
indicated  below. 


0 © J 

3S  + 1 

® 

(input)  1 

.2S3  + 1.2S2  + S 

(output) 

Figure  25.  Unity  Feedback  System 
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I KMNSF  EK  FUl<CTiOl<S 

ul\A  I iCtlLS 

1- 2( 1 ,,U) 

2- 3(3.,  i(  1 • , 0 / .2,3.  1 .2,2,  1 . , 1 I 

3- 2(-  U ,o) 

Note:  All  constant  coefficients  must  include  a decimal  point  as  shown  above. 

Reference  to  any  node  is  made  by  adding  a P to  it.  For  example,  to  obtain 
the  output  at  nodes  2 and  3 the  user  must  refer  to  these  points  as  P2  and 
P3  under  outputs.  A more  complete  input  description  of  Figure  22  will 
clarify.  As  shown  below,  the  input  to  the  system  is  defined  as  PI  and 
described  as  Table  T1 . 


I KANSF Lh  FUNCTIONS 
uKmNCHES 

1-2(1. ,0) 

2- 3(3., 1,  1 . , 0/  .2,3,  1.2,2,  1.,  1) 

3- 2 (- 1 . , 0) 

DEFINED  PARAMETRS 
P1=T  1 
FUNCTIONS 

T1=U,0,  0,1,  10,1 

OUTPUTS 

Pl,P3,PLOT  _ Output  requests  at 

RUN  CUMROLS  Nodes  1 and  3 

STOP  TIME- 10 

INTEGRATION  ROUT INE=IMPLICIT 
END 


The  integer  numbers  between  commas  denote  the  power  of  S associated  with 
the  individual  coefficient. 

The  coefficients  can  be  made  up  of  literal  and  numerical  values. 
Several  sample  branch  transmi ttances  will  better  illustrate  the  required 
i nput : 


Input  at 
Node  1 
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The  branch  descriptions  for  the  flowgraph  shown  in  Figure  Z 3 are  given 
below: 

GRAPH  DESCRIPTION 
BRANCHES 
A-B  (1., 0/1,1) 

B-C  (l.,l/PKl,2,PTAU,l,l,l) 

C-D  (PK2.0) 

C-B  (-PK3.0) 

D-B  (-PK4.0/1.1) 
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SECTION  IV 


DIFFERENTIAL  EQUATION  INPUT  LANGUAGE 

Differential  equations  are  a very  general  form  of  system  definition 
and  all  the  systems  seen  so  far  can  be  described  in  terms  of  these. 

NELSIM  accepts  as  input  sets  of  n*"*1  order  differential  linear  or  nonlinear 
differential  equations  of  the  form  given  below. 


The  above  equations  are  a general  set  of  m nonlinear  differential  equations 
in  m unknowns,  x ^ through  and  with  forcing  functions  F^  through  Fm. 

For  convenience,  the  equations  are  represented  in  matrix  form,  but  will  be 
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input  into  NELSIM  individually. 


The  term 


is  the  nth  derivative 


of  the  variable  raised  to  the  iL  power  where  i can  vary  from  1 to  i 
and  is  not  necessarily  the  same  for  all  derivative  orders.  For  the  systems 
considered  in  Section  II,  the  highest  order  of  derivative  n generally 
encountered  is  two.  The  coefficients  a^  , , etc.  may  be  time  varying  or 

nonlinear.  It  will  be  possible  to  enter  the  coefficients  into  NELSIM  as 
constants,  defined  parameters,  or  any  other  type  documented  in  section. 

NELSIM  automatically  takes  equations  and  separates  them  into  first 
order  differential  equations  in  a format  directly  compatible  with 
SCEPTRE.  It  should  be  realized  that  to  determine  the  transient  response 
of  the  system,  it  will  be  necessary  for  SCEPTRE  to  integrate  a number  of 
equations  equal  to  the  total  number  of  derivatives  - up  to  n derivatives 
for  each  unknown  for  a maximum  of  mn  equations  to  integrate.  The  maxi- 
mum number  of  first  order  differential  equations  currently  allowed  in 
SCEPTRE  is  one  hundred,  and  NELSIM  will  warn  the  user  if  over  one 
hundred  equations  are  generated.  If  the  SCEPTRE  limit  is  increased, 

NELSIM  will  be  easily  modified  to  accommodate  the  increase. 

VARIABLE  PARAMETER  NAMES 

All  names  entered  in  the  descriptions  must  be  entered  as  defined 
parameter  names.  For  example  velocity  VI  in  an  equations  should  be 
expressed  as  PV1 . The  value  of  these  variables  are  then  defined  under 
the  DEFINED  PARAMETERS  subheading. 

SPECIFICATION  OF  DERIVATIVE  ORDER 

The  order  of  a derivative  is  specified  by  a D followed  by  an  integer 

2 

number  indicating  the  derivative  order.  For  example^  d V 2 is  defined  as 

dt7 
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D2PV2  where  PV2  is  the  defined  parameter  name  give  to  V2  and  D2  denotes 
the  second  order  derivative  of  V2. 

As  an  example  of  the  differential  equation  input  in  NELSIM,  consider 
the  equations  below  which  represent  the  two  second  order  differential 
equations 

K1  (V2-V1 ) + M2  V2  + B (V2-V3)  + K2  (V2-V3)  = 0 
B ( V3-V2)  + K2  (V3-V2)  + M3  \/3  + K3  V3  = 0 


resulting  from  the  mechanical  system  shown  in  Figure  10. 

The  NELSIM  input  for  differential  equations  (3)  and  (4)  is  shown  in 
Figure  26. 


U i t h ERENT l ML  EQUATIONS 
UNI  i 3-ENUL ISH  (i  .ECHAiv ICAL ) 

K 1 * ( V2- V 1 ) + M2  * 02  V2  + t * (Q1V2-UIV3)  + K 2 * (V2-V3)  - U 
u * (OIV3-U1V2)  + K2  * C V 3 — V 2 3 + M3  * (u2  V3)  + K3  * V3  - u 
UlT inLU  PARAMETERS 
Kl  = l. 

K2  = 2. 
i\2~. 2 
M3-.3 

U—  1 . 0 

V 1 =X 1 (CGS(2*TIHE)) 

OUTPUTS 
V2 , V3, PLOT 
RUN  CClUlKOLS 

PUhCH  UUIPUT  REQUEST  TO  OBTAIN  DECK  OF  INPUT  CARDS  READY 

STUP  1II.E  = 1U0  FOR  SCEPTRE  INPUT 

ENU 


Figure  26.  NELSIM  Differential  Equations  Input 


73 


The  symbols  D1  and  D2  preceding  a variable  refer  to  the  first  or  second 
time  derivatives  of  that  variable  and  the  defined  parameters  express  the 
variables  involved. 

INPUT  RESTRICTIONS 

To  accommodate  the  SCEPTRE  input  language,  several  restrictions  are 
placed  on  the  input  of  differential  equations  and  are  listed  below: 

1.  The  highest  order  term  must  be 

a)  Uncoupled 

b)  Not  exceed  an  order  of  9 

2.  Derivative  variables  cannot  exceed  4 characters 
(6  including  the  derivative  definition),  e.g., 

D2PV2 

D1PX14 

3.  The  highest  order  derivative  must  appear  only  once  within 
the  equation,  and  be  of  one  of  the  following  forms  only: 

± (expression)  * derivative 
i expression  x derivative 


I 


SECTION  V 


COMBINED  SYSTEM  INPUTS 

A complete  system  is  usually  made  up  of  a number  of  subsystems  not 
necessarily  each  being  of  the  same  type.  NELSIM  accepts  as  input  a 
combination  of  systems  as  well  as  simple  systems  as  seen  in  the  previous 
sections.  A valid  combined  system  entry  of  the  problem  into  NELSIM  is 

ANALOG  DESCRIPTION 

THERMAL 

o 

o 

o 

o 

o 

MECHANICAL 

o 

0 

0 

I 0 

0 

ELECTRO  OPTICAL 
o 
o 
o 
o 
o 

THERMAL 

o 

0 

0 

0 

0 

MECHANICAL 
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SECTION  VI 

NELSIM  SAMPLE  PROBLEMS 


The  following  sample  problems  have  been  formulated  to  illustrate  the 
use  of  NELSIM  with  each  of  the  four  disciplines  considered  - mechanical, 
thermal,  electro-mechanical,  and  electro-optical. 


6.1  MECHANICAL  SYSTEM  EXAMPLES 

Problem  1 - The  problem  depicted  in  Figure  27 


M3  = .3 
B = 1 
V-j  = cos  2t 


Figure  27.  Sample  Mechanical  Problem  1 
can  be  represented  as  shown  in  Figure  28  which  carries  a NELSIM  input 
as  shown  below. 


ANALOG  DESCRIPT  ILL 

MEChANICAL 

ELEMENTS 

VI , L-3-G 1 (TIME) 

K3,G-1=3 

M3, 1 -G=.3 

K2 , 1-2=2 

b 1 ,1-2=1 

M2 , 2-G=. 2 

K1 ,2-3=1 

f UNCTIONS 

Q 1 ( T ) = ( CCS ( 2 . * T ) ) 

OUTPUTS 

VM2(V2),VM3(V3),PLOT 
RUN  CONTROLS 

INTEGRA! ION  ROUT  1 NE  = I MPL I C I T 

£ T Of  T IME=1  0 

END 
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Figure  28. 

Sample  Mechanical  Problem  1 

Representation  to  NELSIM 

The  output  from  NELSIM  will  be: 

CIKCU1 1 

UESCRIP1 IGI 

E UNCT  III* 

S 

QD 1 V ( A)= 

(1./A) 

Q 1 (1 )=(CUS ( 2. *T ) ) 

ELEMENTS 

El 

,0 

-3  =Q! ( T IME ) 

L3 

,c 

-1  = 3. 333330-01 

C3 

,1 

-G  = 3.00000E-01 

L2 

, 1 

-2  = 5.COOOCL-G1 

RI 

, 1 

-2  = 1.C0000E+0L 

C 2 

,2 

-G  = 2 . OOOOOL- G 1 

LI 

,2 

-3  = 1 . OOOOOE+OU 

OUTPUTS 

VC2(V2)f 

PLLT 

VC3(V3), 

PLOT 

RUN  CONTROLS 

STOP  TIME  = 

1 .000E+01 

INTEGRA! I OKRCUT 1 ME 

= IMPLICIT 

Figure  30.  Sample  Mechanical  Problem  2 
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Can  be  represented  to  NELSIM  as: 


ANALOG  LitSCk  IPT  ION 

MECHANICAL 

ELEMENT  S 

W 1 , G- 1 =1 

K 1 ,1-2=1 

J2 , 2-G= 1 

K2 , 2- 3= 1 

J3,3-G=1 

d3  , vj3~G=  1 

L.2,  J2-L-1 

OUTPUTS 

W J2 , UB2 , PLOT 

RUN  CONTROLS 

INTEGRAT ION  ROUT  I NE  = IMPL IC IT 

STOP  T 1/  0 = 1 

END 


The  output  from  NELSIM  will  be 


CIRCUIT  DESCRIPTION 
FUNCT IONS 
QDIV(A)=( 1 ./A) 
ELEMENTS 


E 1 

,G 

-1 

— 

1 .OOOOOE+OC 

LI 

, 1 

-2 

1 . OOOOGE+CU 

C 2 

,2 

-0 

— 

1 .OOOOOE+OC 

L2 

,2 

-3 

— 

1 .OOOOOf  + 0f 

C 3 

,3 

-0 

— 

1 .000000+00 

R3 

,3 

-0 

= 

1 .OOOOOE+OC 

R2 

,2 

-0 

= 

1 .OOOOOL+OC 

OUTPUTS 
VC2.PLCT 
VR2 , PLOT 
RUN  CONTROLS 

STOP  TIME  = 1 . POOE+OO 

INTEGRATI0NR0UT1NE= IMPLICIT 
END 
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6.2  THERMAL  SYSTEM  EXAMPLES 


Problem  1 


Test  Stand 


(a)  Rocket-engine  walls 


(b)  System  Model  Showing 
^ Figure  31.  Thermal  System  Problem  1 
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The  input  of  the  problem  requires  the  various  subdivisions  within  the  wall 
to  be  indicated  by  nodes  (1-3).  The  thermal  resistance  between  subdivisions 
is  simulated  using  resistors  and  thermal  capacitance  using  capacitors. 

The  input  is  shown  below. 


ANALOG  DESCRIPTION 

THERMAL 

ELEMENTS 

TIN ,G-1=7 1 

Rl , 1-2=1 

Cl ,2-0=1 

R2, 2-3=1 

C2, 3-0=1 

R3, 3-4=1 

TLOT ,0-4=1 

FUNCTIONS 

T1 =0,1, I, 1,2, 10, 5, 10, 6, 1,10,1 

GIT PUT  S 

T 1 N , QC 1 , GC2 , QR3 , PLOT 
Ti.3  , TC  1 , PLCiT 
RUN  CONTROLS 

INTEGRA  7 I(  N ROUT  I NE=IMPL I C I T 
STOP  TIME  =15 
i<CN  INITIAL  CONDITIONS 
END 


The  elements  have  all  been  entered  as  constant  values  for  simplicity  but 
in  general  will  be  nonlinear  functions  of  the  materials  involved  and 
depend  on  assumptions  made  by  the  user.  All  the  nonlinear  entry  capa- 
bilities allowed  in  SCEPTRE  apply  to  the  NELSIM  modules.  The  output  of 
NELSIM  is  shown  below  with  the  elements  and  output  requests  altered  to 
electrical  analogies. 
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CIRCUIT  UESCRIP1  I Of. 
FUI.CT  Ii  i.S 
QDIV(A)=( 1 ./A) 

T 1=0, 1,1, 1,2,10,5, 
LL Li. El'll  S 


Ell,  , L 

R 1 , 1 

Cl  ,2 

R2  ,2 

C2  ,3 

R3  ,3 

EOU I ,0 


OUTPUT  E 
Eli.,  PLl  i 
1C  1 .PLOT 
IC2 , PLOT 
IR3 , PLOT 
VR3.PLU 


10,6,1,10,1 


- /. 

-0 

-3 

-L 

-E 


VC  1 , PLOT 
RUi.  CGMKOLS 

STL^  TI/.E  - 1 . 500E+0 1 

I NT  EGRAi  IGNROUT I NE-  IttPL  I C I T 


RUi  INI  I/'  LCOND  I T IONS 
ENL 


= 1 


1 ,000001  4 
1 . 000001  + G( 
1 . OOOOC  E+0< 

1 .OCOQUt  +01 
1 .COOOOL+U 
1 . COOOCE+UO 
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Problem  2 - For  the  rocket  motor  shown  in  Figure  32,  heat  is  transferred 
by  convection  and  radiation  to  the  chamber  wall,  then  through  the  wall  by 
conduction,  and  finally  to  the  coolant  by  convection 


A-A 


Figure  32.  Thermal  System  Problem  2 
This  problem  is  described  to  NELSIM  as  follows: 


ANALOG  DESCRIPTION 

1 HEKP.AL 

ELEMENTS 

7 GAS , G — 1 — T 1 

R1R, 1-2=1 

R1C, 1-2=1 

R2K, 2-3= 1 

C2 , 3-G= 1 

R3C,3-A=1 

TCOGL,G-l*=1 

FUNCT  IGI.S 

1 1 =L , 1,1,1 ,2 , 10, S, 10,6, 1,10,1 
OUTPUT 

T GAS , TCLLL , I LOT 
TR2K, PLOT 
RUN  CONTROLS 

INTEGRAT ION  ROUT  I NE  = IMPLICIT 

STOP  T II  £=1  0 

END 


The  produced  NELSIM  output  for  this  case  will  be; 
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CIRCUIT  DESCR  IPT  IGf 
ELi.CT  II  . S 
C D I V ( A ) - ( 1 ./A) 


T 1 =t .1.1,1 

,2, 

10,5, 10,6, 1,10,1 

LLLI  ENTS 

EGAS 

» G 

- 1 

=T  1 

K 1 1 v 

, 1 

1 . 00000E+  GO 

RIC 

, 1 

r 

— 

] .QOOGOi  +0( 

R2K 

,2 

-3 

— 

1 . OOOOCL+OG 

C 2 

,3 

-G 

— 

1 . 000001  +0C 

R3L 

,3 

-i* 

1 . OOOOOL+CG 

ECOCL 

» G 

-k 

- 

1 . 000001  +0G 

OUTPUTS 

EGAS, PLOT 

ECGOL , 1 LOT 

VR2K, PLOT 

RUN  CGNTROLS 

STLP  T li  t. 

= 1.C0CL+01 

1 i.TL  GK/  T I Ci  .ROUT  I RE  - 1 MPL  I C I T 
END 
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6.3  ELECTRO-MECHANICAL  EXAMPLES 

Problem  1 - For  the  galvanometer  shown  in  Figure  33. 


Figure  33.  Electro-Mechanical  Problem  1 

Figure  34  represents  the  equivalent  lumped  network  in  terms  of  its  elec- 
trical and  mechanical  parts. 


Figure  34.  Lumped  Element  Galvanometer  Equivalent 
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where  J and  B represent  the  moment  of  inertia  of  the  rotor  and  the  viscous 
damping  that  results  from  air  friction.  The  element  represents  the 
rotational  stiffness  produced  by  a spring  attached  to  the  rotor.  The 
constant  D is  defined  to  be  eta  where  8 is  the  magnetic  flux  density,  i is 
the  total  length  of  the  conductor  and  a is  the  radius  of  the  rotor. 

The  input  language  where  all  the  elements  are  normalized  to  one  is 
shown  below. 


ANALOG  DESCRIPTION 
ELECTRO  MECHANICAL 
ELEMENT S ( I LECTRICAL) 

El , G- 1 - 1 
R 1 , 1-2-1 

EC,G-2=CEC(PD,WJR) 

ELEMENT  S (MECHANICAL ) 

T 1 , G-3-CLC ( PO , I R 1 ) 

JP. , 3-G=  1 
BE , 3-0=1 
KRf3-G=! 

DEFINED  PARAMETERS 
PD=QPO(PA,PL,PB) 

PA=1 
PL—  1 
PC=1 
PL=  1 

FLNCT  It  NS 
QEC(A,u)=(A*8) 

QPD ( A, u , C ) =( A*B*C ) 

GLTPUTS 

t.KR(KKR), TKR(TKR), VK1 .PLOT 
RON  CGiH  HOLS 
STOP  T 1 ME  = 1 0 

INTEGRATION  ROOT  I NE  = IMPL I C 1 7 
END 


The  two  dependent  sources  EC  and  T1  provide  the  interface  between  the  two 
systems  and  are  entered  directly  in  terms  of  the  elements  involved.  The 
output  of  NELSIM  is  shown  below  with  the  elect* ical  analog  represented 
shown  in  Figure  35.  As  can  be  seen,  the  network  is  entirely  in  terms  of 
electrical  quantities. 
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CikCUn  UtbCR  IP  i I ui. 
F Ui.CT  I ( AS 


GD1V(A)= 

( 1 ./A) 

QEC(A.L) 

=(/ *B) 

QPL  (A,  t. , 

C ) = ( A*b';  C ) 

ELEMENTS 

J1 

,C  -3 

=QLC(PD, IR1 ) 

CR 

,3  -C 

= 1 . OPCODE + 0; 

KR 

,3  -C 

= 1.0000CL  + D- 

LR 

,3  -0 

= 1.COOOCMC 

E 1 

-1 

= 1 . ODOOOL+b 

R 1 

, 1 -2 

1 . COCOOL+O 

EC 

,(■  -2 

=QEC(PL, VCR) 

DEE  IREL 

PARAMETERS 

PD 

-CPD ( E A, f L , PD ) 

PA 

= 1 .OOCOOE+OO 

Pb 

I . COCCOE+OO 

PC 

1 . OCiLOCE  + OO 

PL 

= 1 . OOOOOE+OO 

OUTPUTS 
VLR(bj(R),PLOT 
I LF.  ( T KR  ) , PLOT 
VR 1 , PLLT 
RUI.  COM ROLS 

ST LP  T 1/it  = 1 . OOCjE  + O 1 

INI  LGRAT  IU.RCUT  INt=  IF.PL  ICI T 
ERD 


Problem  2 - Implementation  of  the  accelerometer  model  built-in  to  the 
electro-mechanical  system  is  accomplished  below.  The  usage  and  response 
of  the  model  is  illustrated  for  the  case  where  all  gains  have  been 
normalized  to  1 and  the  input  to  the  model  is  a unit  step.  Consider 
Figure  below  where  H is  the  accelerometer. 


© 


The  input  to  NELSIM  is  given  below: 


ANALOG  DESCRIPTION 
ELECTfi-  ML  CHAM  CAL 
MODELS 

T 1 , I N-i  L=/  CCELf  -REMOTER 
OLE  I NEC  P/  KAMETERS 

pa; i=i 

PIT  1 = 1 

PCT 1=1 

PDT 1=1 

PET  1 = 1 

Pf  T 1 = 1 

PoT 1 = 1 

PI  NT  1=1  1 

Fl'NCT  IONS 

T 1=0, 0,0, 1,10,1 

01 T PUTS 

P INTI, PA  I i 1 , PA2T 1 , PA3T 1 , PAA T 1 , PLOT 
PAST  1 , POL i 1 .PLOT 
RLN  COM Rl  LS 

INTEGRA! ION  ROUT  I NE  = IMPL IC I T 

S10P  TIME- 10 

EI.D 


All  gains  within  the  model  are  provided  by  the  user  wtf 
(Tl)  added  to  each  name. 
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' AD-A033  606 


UNCLASSIFIEO 


TRW  SYSTEMS  6R0UP  REOONDO  BEACH  CALIF  NUCLEAR  SURVIV— ETC  F/B  9/2 
NONELECTRICAL  LAN6UAGFS  SIMULATION  MODULE  (NELSIM).  VOLUME  II, (U) 
SEP  76  M EPSTEIN*  J R PISTACCHI  F29601-73-C-0024 


AFWL-TR-73-I56-V0L-* 


NL 


2of2 

10 

A033606 


END 

DATE 

FILMED 

2—77 


V 


The  output  generated  by  the  program  is  shown  below. 


CIRCUIT  DESCRIPTOR 

FURCT  IRRS 

CD  I V ( A)  = ( 1 ./A) 

T 1 =0, 0 , 0 , 1 ,10,1 

QMaI (A,D,C)=(A*l-C) 

QMA2 ( A, b , C , D )=( ( A*L-C ) /D ) 

QMA3 ( A,b )=( A*B) 

OM.A  L ( A , c , C ) = ( ( A- u ) / C ) 

CM/  S(A,B,C,D)=(>  * ( ( L.-C  ) /D  ) ) 

CM/ 6( A , L ) =( A/B ) 

DETiliEL  PARAMETERS 
PA1  1 = 1 . OOC.  OOE+OO 

PUT  1 - l.OClCOE+OO 

PCT 1 - 1.00( OOE+OO 

PDT 1 = 1. OOC OOE+OO 

PET  1 = 1. COT OOE+OO 

PET  1 = l.OOCOOE+OO 

P011  = 1 . GOCOOE+OO 

PIRT1  = T 1 

PA  1 T 1 =C,MA  1 ( P I NT  1 , P/  T 1 , PAA  T 1 ) 

QP/.2T  !=PA7T  1 

DPA7T  1=CMiA2  (PA IT  1 ,f  fcT  1 , PA7T  1 , PCT  1 ) 
PA3T 1=CMA3 (PA2T 1 , PL  T 1 ) 

DPA6T  1 =GAiAA  ( PA3  • 1 , i A6T  1 , PP  T 1 ) 

PA*»T  1=CMA5(PET  1 ,PA3T  1 ,PA6T1  ,PFT  1 ) 
DPAST  1 =CNA6( PALI  1 ,P(.T  1 ) 

DPLUT 1=PA5T 1 
PA2T 1=0 
PAST  1=0 
PA6T 1=0 
PA7T 1=0 
POUT  1=0 
OUTPUTS 
PIRT 1 .PLOT 
PA  1 T 1 .PLOT 
PA2T1 .PLOT 
PA3T 1 , PLOT 
PALT1 .PLOT 
PAST  1 .PLOT 
PCLT1 .PLOT 
RUR  CORTKOLS 

STLP  TIME  = 1 . UOOE+O 1 

I NTEORAT ICRROUT 1RE=1MPL IC I T 
END 
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6.4  ELECTRO-OPTICAL  SYSTEM  EXAMPLES 

Figure  36  details  the  equivalent  circuit  selected  to  simulate  the 
AC  characteristics  of  photo  diodes. 


/ 


Figure  36.  Equivalent  Circuit  of  a Photo  Diode 


As  detailed  In  Volume  I,  the  photocurrent,  IpH  is  expressed 
as 

IpH  • P/hV) 

Problem  1 

The  normalized  power  input  at  5 m to  an  INSb  device  is  described  as 
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Using 


q = 25% 

Rs  = 18  n 
Cf  = 7.1  pf 
R ^ s S FI 

rload  = 12 n 

Plot  the  current  time  history  at  the  load  and  at  the  source. 
The  NELSIM  input  is 


ANALOG  DESCRIPTION 
ELECTRO  OPTICAL 
ELEMENTS (ELECTRICAL ) 

JPD , G- 1=CEO(PEFF , PChG, PCON, POFR, TPL ) 

R I , 2-G=6 

C 1 , 3-G=3 

RS, 3-4=12 

RL , 4-G= 1 2 

FUNCT IONS 

TPL=C,0,.2,2,.4,6,.E5, 10,.  6,  12,0,  15 

CEO(A,b,C,D,E)=(A*B/(C*D)*E) 

DEFINED  PARAMETERS 

PEFE=.2S 

PCH(.=  1 E-7 

PCON= 1 23E-4 

POFP=37 

COT RUTS 

JPD, JRL, PLOT 

RUN  CONTROLS 

INTEGRA! ION  R0UT1NE=IMPLICIT 
STOP  T IME=1  S> 

END 
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For  which  the  NELSIM  output  Is 


CIRCUIT  DESCRIPT  ION 

FUNCT  li  NS 

CD  I V ( A ) = ( 1 ./A) 

TPL=0,O,.2,2,.4,6f .85, 10,. 6, 12,0, 15 
QEL ( A , b , C , D , E )=( A*L / ( C*0 ) *L ) 
ELEMENTS 

JPC  ,G  -1  =QE 

RI  ,2  -0  = 

Cl  ,3  -C  = 

RS  ,3  -A  = 

RL  , ^ -G  = 

DEFINED  PARAMETERS 


PEFF  = 

PCFiG  = 

PCCN  = 

PCF  R = 

OUTPUTS 
JPC, PLOT 
JPxL , PLOT 
RUt.  CONTROLS 
STOP  TII  E 


2.50C00E-0I 
1 . C0000E-07 
1 .23C00E-02 
3.7Q(.OOE+Ol 


=QEO( PEFF ,PCHC, PCOK, Pf  f K# Tf L ) 
= 6. POOOOL+Ot 

= 3.  C 00  ODE  + 0( 

= 1.20000E+01 

= 1 . 20000L+D 1 


1.500E+01 


I NT  ECRAT  IOf  .ROUT  INE=  IMPLICIT 
END 


Problem  2 


Rerun  Problem  T using  an  INA$  device  whose  characteristics  are 


q = 25% 
Rs  = T2  n 


For  the  same  Toad  and  compare  results. 
The  NELSIM  input  and  output  are 
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ANALOG  DESCRIPTION 
ELECTRO  OPTICAL 
ELEMENTS (ELECTRICAL) 

JPD,G-  1=GEG(PEFF,  PCliG,  PCON,F'CFR,  TPL ) 
R 1 , 2-G=6 
C 1 , 3-G=7. 1 
RS, 3-4=10 
RL , 4-G=1 2 
t UNCT ICRS 

TPL=G,0,.2,2,.4,6,.E5, IV, .6, 12,0, 15 
CEO(A,B,C,D,E)=(A*D/(C*D)*E) 

DEFINED  PARAMETERS 

PEFF=. 25 

PCHC= 1 E- 7 

PCON=123L-4 

POFR=37 

OUTPUTS 

JPD, JRL , PLOT 

RUN  CONTROLS 

INTEGRA! ION  ROUT  1 NE  = IMPL IC I T 

STOP  T I ME = 1 5 

END 


CIRCUIT  DESCRIPT  ION 
FUNCTIONS 


QD I V ( A)= 

(1./A) 

TPL=0,G, 

.2,2, .4,6,. 85,10,. 6, 

12,0,15 

QEL( A,b, 

C,D,E)=(A*L/(C*D)*E) 

ELEMENTS 

JPD 

,G  -1 

=QEO(PEFF,PCHC, 

R I 

,2  -G 

= 6. OOOOGE+OG 

Cl 

,3  -G 

= 7.  1 OOOOE+OO 

RS 

.3  -4 

= 1 . EOOOOE+0 1 

RL 

,4  -G 

= 1 . 2000GL+G1 

DEFINED 

PARAMETERS 

PEFF 

= 2. 50000E-01 

PCPiG 

= 1 . 00G00E-07 

PCCN 

= 1.23C00E-02 

POFR 

= 3. 70C00E+0 1 

OUTPUTS 
JPD , PLOT 
JRL, PLOT 
RUN  CONTROLS 
STOP  TIME 


1 . 500E+0 1 


INTEGRAT IONROUT IRE=IMPLICIT 
END 
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